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Abstract. We survey results that have been obtained for self-adjoint op- 
erators, and especially Schrodinger operators, associated with mathematical 
models of quasicrystals. After presenting general results that hold in arbitrary 
dimensions, we focus our attention on the one-dimensional case, and in par- 
ticular on several key examples. The most prominent of these is the Fibonacci 
Hamiltonian, for which much is known by now and to which an entire section 
is devoted here. Other examples that are discussed in detail are given by the 
more general class of Schrodinger operators with Sturmian potentials. We put 
some emphasis on the methods that have been introduced quite recently in the 
study of these operators, many of them coming from hyperbolic dynamics. We 
conclude with a multitude of numerical calculations that illustrate the validity 
of the known rigorous results and suggest conjectures for further exploration. 
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1. Introduction 



The area of mathematical quasicrystals is fascinating due to its richness and 
very broad scope. One perspective that leads to a rich theory is the question of 
how well quantum wave packets can travel in a quasicrystalline medium. This, in 
turn, leads to a study of the time-dependent Schrodinger equation governed by a 
potential that reflects the aperiodic order of the environment to which the quantum 
state is exposed. 

This survey paper describes the current "state of the art" of mathematical re- 
sults concerning quantum transport in mathematical quasicrystal models. We will 
describe the (Schrodinger) operators that are typically considered in this context, as 
well as the spectral and quantum dynamical properties of these operators that are 
relevant to our basic motivating question. Few results hold without further assump- 
tions. After presenting these general results, we will therefore specialize our class 
of operators in various ways. On the one hand, additional general results follow 
from passing to a one-dimensional setting, most notably some specific consequences 
of Kotani theory. On the other hand, there is much interest in certain central ex- 
amples, such as the Fibonacci model, and more generally the Sturmian case or 
potentials generated by primitive substitutions. For these special cases, many more 
results are known, most of which will be described in detail here. Another purpose 
of this survey is to highlight new tools that have recently been introduced in the 
study of these models that have led to very fine quantitative results in the Fibonacci 
case and beyond. 

The structure of the paper is as follows. In Section [2] we describe the operators 
that are typically considered in the case of general dimension. In Section [3] we 
present the results known to hold in this very general setting. Next, in Section^ we 
pass to the one-dimensional situation, where the operators are slightly redefined to 
conform with the bulk of the literature. Additional general results are described in 
this scenario. These range from Kotani theory, via proofs of zero-measure spectrum 
and hence absence of absolutely continuous spectrum based on Kotani theory, to 
proofs of absence of point spectrum and hence purely singular continuous spectrum 
based on Gordon's lemma. We also discuss how these general results can be applied 
to several classes of examples. Section [5] is devoted to the special case of the 
Fibonacci Hamiltonian. This is the most prominent one-dimensional quasicrystal 
model; in addition to the general results mentioned before, there are fine estimates 



SCHRODINGER OPERATORS ARISING IN THE STUDY OF QUASICRYSTALS 3 

of the local and global dimensions of the spectrum and the density of states measure, 
as well as the optimal Holder exponent of the integrated density of states, at least in 
the regimes of small or large coupling. For this model the question about quantum 
transport behavior has very satisfactory answers, at least in the regime of large 
coupling. All these results, along with comments on their proofs, are presented in 
this section. Then, in Section |6j we discuss how the approach and the results may 
be extended from the Fibonacci case to the more general Sturmian case. We present 
and discuss numerical results in Section [7] Finally, Section |8] lists and discusses a 
number of open problems that are suggested by the existing results. 

Related matters have been surveyed earlier in (TUl SU SH 1147] . 

2. SCHRODINGER OPERATORS ARISING IN THE StUDY OF QuASICRYSTALS 

In this section we describe self-adjoint operators that have been studied in the 
context of quasicrystals. There is a clear distinction between the case of one space 
dimension and the case of higher space dimensions. While the geometry of the 
quasicrystal model plays an important role in higher dimensions, this is not the 
case for one-dimensional models. In the former case, this leads to a dependence of 
the underlying Hilbert space on the realization of the model. In particular, as one 
typically embeds any such realization in a family of realizations, this leads to some 
technical issues that need to be addressed mathematically. In the latter case, on 
the other hand, one usually works in a universal Hilbert space and the aperiodic 
order features are solely reflected by the potential of the operator. This allows one 
to invoke the standard theory of ergodic Schrodinger operators with a fixed Hilbert 
space. In this survey we will present the known results in the settings in which they 
have been obtained, which we now describe. 

Quasicrystals are commonly modeled either by Delone sets in Euclidean space 
or by tilings of Euclidean space. In fact, any such Delone set or tiling is embedded 
in a Delone or tiling dynamical system, which is obtained by considering the set 
of translates and then taking the closure of this set with respect to a suitable 
topology. This orbit closure is called the hull of the initial Delone set or tiling. 
The dynamics is then given by the natural action of the Euclidean space on the 
hull through translations. For this dynamical system, one then identifies ergodic 
measures, and they are typically unique. While each element of the hull gives 
rise to a Schrodinger operator, it is the ergodic framework that allows one to prove 
statements that hold for almost all such operators, as opposed to results for a single 
such operator. Occasionally, it is then even possible to extend results that hold for 
almost all elements of the hull to all elements of the hull by approximation. 

Delone sets and tilings are in some sense dual and hence equivalent to each 
otherj^ For definiteness, we will consider a framework based on Delone sets. Let 
us fix a dimension d > 1. The Euclidean norm on M'' will be denoted by | • |. We 
denote by B{x,r) the closed ball in R'^ that is centered at x and has radius r. 

Definition 2.1. A set A C M"^ is called a Delone set if there are r, R > such that 

B{x,r)r\A^ {x} Vx e A 

^One can go back and forth between these two settings by decorations and the Voronoi 
construction. 
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and 

B(a;,i?) n A 7^ VxeM'*. 

Thus for a Delone set we have a lower bound on the separation between points 
of the set, and an upper bounds for the size of "holes" in the set. One also says 
that a Delone set is uniformly discrete and relatively dense. Since Delone sets are 
closed and we want to take orbit closures, let us now define the underlying topology 
on F{W^), the closed subsets of W^, following [TTT]. 

Definition 2.2. For k e Z+ and F,G E we define 

dk{F, G) := inf ({e > : F n B{0, k) C U^SG) and G n 5(0, k) C U^{F)} U {1}) . 

For k £ Z+, e > 0, and F e J'(M''), we define 

U,^k{F) := {G e TiR") : dk{F, G) < e}. 

The topology on J-{M.'^) with neighborhood basis {U^,k{F)} will be called the natural 
topology and denoted by Tnat • 

Proposition 2.3. (a) Translations are continuous with respect to Tnat- 

(b) J-'(M'') endowed with Tnat is compact. 

(c) Tnat is metrizable. 

See for a metric that induces Tnat- 

Definition 2.4. Let n be a set of Delone sets and denote by T the translation 
action ofM.'^, that is, TfX = x + t. The pair (f2,T) is a Delone dynamical system 
if CI is invariant under T and closed in the natural topology. 

Definition 2.5. A Delone dynamical system is said to be of finite local complexity 
if for every radius s > 0, there is a uniform upper bound on the number of different 
patterns one can observe in uj intersected with a ball of radius s. Here, lj ranges 
over n, the center of the ball ranges over , and patterns are understood modulo 
translations. 

Definition 2.6. Suppose {^,T) is a Delone dynamical system of finite local com- 
plexity. A family {A^jj^gji of bounded operators A^ : i'^{uj) — ^^(w) is said to 
have finite range if there is s > such that for all uj £ CI, Ai^{x,y) only depends 
on the pattern of uj in the s -neighborhood of x and y, and Ai^{x,y) — whenever 
\x-y\> s. 

The class of operators so defined encompasses all discrete operators that are 
usually considered as quantum Hamiltonians in the context of multi-dimensional 
quasicrystals. In one dimension, however, it is customary to realign points as a 
lattice (i.e., Z) and to encode the geometry in the matrix elements of the operator. 
Even more specifically, one focuses on nearest neighbor interactions and hence ob- 
tains a tridiagonal matrix in the standard basis of £^(Z). Much of the mathematical 
literature focuses on the discrete Schrodinger case, where the terms on the first off- 
diagonals are all equal to one and the quasicrystalline structure of the environment 
is reflected in the terms on the diagonal. That is, the family of discrete Schrodinger 
operators one then considers is of the form {i?tj}a;Gf2j where f2 is typically a subshift 
over a finite alphabet, and for w £ CI, acts on vectors from £^(Z) as follows, 

[H^^l^]{n) = ij{n + 1) + i,{n - 1) + V^{n)^{n). 
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Here the potential VL, is given by K;(n) = /(r"w), where T : 51 — > is the standard 
shift transformation and / is at least continuous, usually locally constant, and often 
just depends on a single entry of the sequence w. 

3. General Results in Arbitrary Dimension 

3.1. Spectrum and Spectral Types. One of the fundamental results for the 
families of operators introduced above is the almost sure constancy of the spectrum 
and the spectral type with respect to an ergodic measure p associated with (fi, T). 
This follows from the covariance condition 

where Ut : ^^(w) — > £'^{uj + t) is the unitary operator induced by translation by i, 
along with the definition of ergodicity applied to traces of spectral projections in 
the usual way. This establishes the following result; compare [lllj . 

Theorem 3.1. Suppose (O, T) is a Delone dynamical system of finite local complex- 
ity and fjL is an ergodic measure. Let {^c^jwesi be a family of bounded self-adjoint 
operators of finite range. Then there exist S, Spp, Sscj ^ac '"^'^ o, subset f2o C f2 o/ 
full ii-measure such that for every u e ilo; we have ai^A^^) = E, crp.p{A^) ~ Spp, 

Csc(^w) = 5]sc, and (Tac(^[^) = 5]ac. 

3.2. Existence of the Integrated Density of States. Suppose {A^jj^go is a 
family of bounded self-adjoint operators of finite range. For Q C M'^ bounded, the 
restriction A^\q defined on t'^{Q n a;) has finite rank. Therefore, 

n{A^,Q){E) :— ^{eigenvalues of A^\Q that are < E} 

is finite and E i— > ^^ji{Ai^^Q){E) is the distribution function of the measure Pq" 
defined by 

j ^dp^- - _l_Tr(^(^^|Q)), ^ e a(R). 

For s > and Q C M"*, denote by dsQ the set of points in whose distance 
from the boundary of Q is bounded by s. A sequence {Qk\ of bounded subsets of 
is called a van Hove sequence if ^°^a\{Qk) ~^ as fc — >■ cx) for every s > 0. 
The following result was shown in [112) (compare also the earlier papers [571I5H] ). 

Theorem 3.2. Suppose (fi, T) is a strictly ergodic Delone dynamical system of 
finite local complexity. Let {^cjjwen be a family of hounded self-adjoint operators of 
finite range and let {Qk\ be a van Hove sequence. Then, ask ^ oo, the distributions 
of E 1-^ Pq'^{{—oo, E]) converge to the distribution E i— p"^{{—oo, E]) with respect 
to II • lloo; o,nd the convergence is uniform m cj G fl. 

In fact, the limiting distribution can be given in closed form; see |112j . It is 
called the integrated density of states, and the associated measure is called the 
density of states measure. The remarkable feature of this result is the strength 
of the convergence, in that the distribution functions converge uniformly in the 
II ■ lloo topology. This is of particular interest in cases when the limiting distribution 
function has jumps. The next subsection shows that the latter phenomenon may 
actually happen. 
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3.3. Locally Supported Eigenfunctions and Discontinuities of the IDS. In 

dimensions strictly greater than one, the local structure of a Delone set may be cho- 
sen such that suitable finite range operators have finitely supported eigenfunctions 
at a suitable energy. If these local configurations occur sufficiently regularly, it 
follows that the energy in question will be a point of discontinuity of the integrated 
density of states. This observation may now be supplemented in two ways. On the 
one hand, a given Delone dynamical system may be transformed into one that is 
equivalent to the original one in the sense of mutual local derivability, which does 
have the required local configurations. On the other hand, any discontinuity of 
the integrated density of states must arise in this way, that is, through the regular 
occurrence of finitely supported eigenfunctions at the energy in question. These 
issues were discussed in the paper [102] . Let us state the results from that paper 
precisely. 

Theorem 3.3. Suppose (f2, T) is a strictly ergodic Delone dynamical system of 
finite local complexity. Let {Ai^}i^^q be a family of bounded self-adjoint operators 
of finite range. Then there exist (f2',T) and {A'^}i^^fii such that (r2,T) and (51', T) 
are mutually locally derivable and A'^ has locally supported eigenfunctions with the 
same eigenvalue for every a; G 51'. Moreover, A'^ can be chosen to be the nearest 
neighbor Laplacian of a suitable graph. 

Theorem 3.4. Suppose (f2, T) is a strictly ergodic Delone dynamical system of 
finite local complexity. Let {Ai^}^^q be a family of bounded self-adjoint operators 
of finite range. Then E £ R is a point of discontinuity of if and only if there 
exists a locally supported eigenfunction of A^ to E for one {equivalently, all) G il. 

4. General Results in One Dimension 

Starting with this section, we will focus on the case of one space dimension. Far 
more rigorous results are known for this special case than for the general case. In 
particular, much is known about the structure of the spectrum as a set, as well as 
the type of the spectral measures. 

As mentioned above, in the one dimensional setting one typically passes to a 
somewhat different choice of the model. Thus, for definiteness, we will restrict our 
attention in much of the remainder of this paper to the following scenario. We 
consider Schrodingcr operators in 

(1) [H^i^]{n) = i:{n + 1) + ^j{n - 1) + K,WV^(n), 

where the potentials are of the form f{T"'uj), with a; in a compact metric 

space il, a homeomorphism T : fl — > Jl, and / G C(il,M). We also fix an ergodic 
measure /u. 

Notice that the Hilbert space in which Hi^ acts is now w- independent, and the 
aperiodic order features of the medium that is being modeled are completely sub- 
sumed in the potential of the operator H^^. In this we follow the standard 
convention, for this class of operators has been commonly studied. One could con- 
sider operators that are formally more akin to the operators considered above in 
general dimensions. However, this would not lead to any significant mathematical 
difference. Loosely speaking, the aperiodically-ordered Delone set in R is just be- 
ing reconfigured as Z, and the local properties of the operator that depend on the 
pattern near a point in the general setting affect the value of the potential at the 
point in question accordingly in our present setting. 
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In fact, this scenario is more general than considered in the previous section. 
The better analog would be the case where {Tuj)„ — con+i is the shift on for 
some finite set ^, C A^ is T-invariant and closed, and / : — M is locally 
constant, that is, it only depends on the values of ujn for some finite set of n values. 
However, some of the results below hold in the more general setting, and we will 
impose further restrictions when they are needed. 

Consequently, Theorems |3.1| and |3.2| now take the following form. 



Theorem 4.1. There are sets S, Sac, "^sc, o-nd Spp such that for fi-almost every 
uj eil, we have a{H^) = S, crac(^^w) = ^ac, f^sc{H^) = Egc, and (Jpp{H^) = Epp. 

Theorem 4.2. The measures j ipdN'^ — ^TT{(f[H^\^i^^)) converge weakly to the 
measure J ipdN = J{So,(p{H^)Sq) dfj,{uj) as k oo. 



The second result uses a weaker notion of convergence than in Theorem |3.2[ the 
price we have to pay for casting this problem in the more general setting. However, 
this is fine after all, due to the following result. 

Theorem 4.3. The measure dN is continuous. 

Moreover, we have: 

Theorem 4.4. The topological support of the measure dN is equal to E. 



Theorems |4.1f[4.4| are standard results from the theory of ergodic Schrodinger 
operators on ^^(Z); compare [50] . 

4.1. Spectrum and the Absence of Uniform Hyperbolicity. Let us consider 
solutions to the difference equation 

(2) u{n + l) + u{n- 1) + V^{n)u{n) = Eu{n), n e Z, 

for some energy E Clearly, u solves ^ if and only if 

Since — f{T^uj), one naturally defines 

SO that ^ implies 



E-fiio) -1 
1 



uin) ')=AEiT-^)-----AE{Tu) (^^J^j 

for n>l and solutions u to ([2]). We set ME,uj{n) = ^£;(r"w) x • • • x ^^(Tw). 
Definition 4.5. We let 

Un = {E ^M.:3c>l such that\/uj e f7, n > 1 we have \\ME,uj{n)\\ > c"}. 

Johnson showed in [57] that the set UH is equal to the resolvent set of for 
any lu that has a dense T-orbit. As a particular consequence, we may state the 
following: 

Theorem 4.6. Suppose T is minimal. Then for every uj € Q, we have a{Hi^) — 
M \ WH. In particular, for any ergodic measure fj,, we have S = M \ WH. 
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4.2. Kotani Theory. In the previous subsection, we saw that the partition M = 
UHLiR\UH yields the partition of the energy axis into resolvent set and spectrum. 
Let us partition the energy axis even further by introducing the Lyapunov exponent 

L^iE) = lim - / log\\ME,Un)\\ d^i{uj). 

n—foo n J 

The existence of the limit follows quickly by subadditivity. Moreover, due to er- 
godicity of fi and Kingman's Subadditive Ergodic Theorem, we have 

Ln{E) = lim — log II Af^; for /i-a.e. uj & ft. 

n— foo 77, ' 

Obviously, we have Lf^{E) > for every E e UH. We let 

Z^ = {E:L,,iE)=0} 

and 

AfUH^, = {EeR: L^,{E) > 0}\UH, 
so that our final partition is M = UT-LuMWH^UZ^. Notice that in these definitions, 
the dependence of the partition of M \ WH into J\fUHf_i U on the choice of the 
ergodic measure ^ is emphasized through the subscript for the /^-dependent sets. 
It is customary to drop this explicit subscript from L, Z, and NWH, and we will 
henceforth do so as well. 

Recall that the essential closure of a set 5 C M is given by 

5°"" = {i; e M : Ve > we have Leb(5 {E - e, E + e)) > 0}. 

Then, the following result is the celebrated Ishii-Pastur-Kotani theorem; see, for 
example, [151 [Ml [TOm im [Ti^ . 

Theorem 4.7. We have Eac = 'Z^'''' ■ 

Since the essential closure of a set S is empty if and only S has zero Lebesgue 
measure, this result yields a characterization of almost sure purely singular spec- 
trum. In fact, this is the typical situation for our models of interest |107| . 

Theorem 4.8. Suppose the range of f : fl ^ R is finite and Leb(Z) > 0. Denote 
the push-forward of fj, under fl 3 lu i-^ Vui G (Ran /)^ by fi* . Then, supp /i* is 
finite. 

Since supp /i is T-invariant, this means that every element of supp /i* is periodic. 
In other words, if the potentials are ergodic, aperiodic and take finitely many val- 
ues, then Z has zero Lebesgue measure and the almost sure absolutely continuous 
spectrum is empty. 

4.3. Zero-Measure Spectrum. The realization that Leb(-Z) = for potentials 
that are ergodic, aperiodic and take finitely many values is at the heart of proofs 



of zero-measure spectrum in these cases. Whenever Theorem 4.8 applies, it suffices 
to show that the spectrum is contained in (and hence coincides with) Z in order 
to establish zero-measure spectrum. There are two approaches that establish this 
identity. One shows that the Lyapunov exponent must vanish for every energy in 
the spectrum. The other is based on Theorem |4.6( which says that E = ZUAfUH. 
Hence, it suffices to prove that if E is such that the Lyapunov exponent is positive, 
the convergence of ^ log ||M£;_^(n)|| to L{E) must be uniform. The second approach 
can be made to work in greater generality, whereas the first approach often gives 
additional information that has other interesting consequences. 
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We will discuss how the first approach is implemented when we discuss Sturmian 
models and, specifically, the Fibonacci Hamiltonian in later sections. Here we 
therefore discuss in some detail how the second approach works. We first state the 
main result, then explain how it is a natural consequence of the line of reasoning 
employed. We emphasize that this is a result that holds in the symbolic setting. 
That is, for a finite set A, called the alphabet, we consider the shift transformation 
T : given by {Tuj){n) = a;(n+ 1) and closed, T-invariant subsets ^ of A?', 

which are called suhshifts [over A). We denote by Wji the set of all finite words 
over the alphabet A that occur in elements of il, and we write Wn{n) for those 
elements of Wn that have length n. A function on 51 is called locally constant if 
f{(jj) only depends on finitely many entries. More formally, there exists k G Z+ 
such that /(w) = /(w') for all w, w' G 51 with ui_k . . .uj^ = ui'_^ ■ ■ ■'^'k- 

As usual, a subshift 51 is called minimal if the topological dynamical system 
(51, T) is minimal (i.e., the T-orbit of every w G 51 is dense in 51), it is called 
uniquely ergodic if (51, T) has a unique invariant Borel probability measure, and it 
is called strictly ergodic if it is both minimal and uniquely ergodic. In the uniquely 
ergodic case, the unique invariant Borel probability measure must necessarily be 
ergodic. 

Definition 4.9. Let be a strictly ergodic subshift with unique T-invariant mea- 
sure /X. It satisfies the Boshernitzan condition (B) if 

limsup min n ■ ^ {[w]) ) > 0. 

Here is the result from [SB] showing that (B) is a sufficient condition for zero- 
measure spectrum: 

Theorem 4.10. Suppose VL is a strictly ergodic subshift that satisfies condition (B) 
and f : Vl ^ M. is locally constant. Then the convergence of ^\og\\M e ^^oin)]] to 
L{E) is uniform for every i5 G M. In particular, AfWH is empty and Yj — Z . Thus, 
if fl and f are such that the T4; are aperiodic, then Leb(S) = 0. In the latter case, 

Sac = 0. 

The proof of Theorem |4.10| is more easily understood if one imposes a stronger 
assumption. Indeed, let us assume for the time being that /(w) depends only on ujq 
and minj„gvVj^(„) {[w]) is uniformly large for all n, not merely for a subsequence. 
That is, suppose there is (5 > such that \w\ ■ ^ {[w]) > S for every w G Wo- 

By our stronger assumption on /, we may view log ||M£;^(^(n)|| as a quantity 
associated with the word w = Wi . . .a;„ G yVn(n) that we will denote by F{w). If 
we do so, then the goal is to prove that \w\~^ F(w) converges uniformly as \w\ oo 
and each w belongs to Wn- It is a known consequence of unique ergodicity of (51, T) 
that the convergence in 

F'^ :— limsup \w\~^F{w) 

is uniform. By the uniform upper bound and the frequent occurrence of any word 
due to the strengthened assumption, one can derive a similar uniform result for the 
liminf, and hence establish uniform convergence. If one only has condition (B), 
then a similar way of reasoning can be carried out for the sequence of length scales 
from (B). In this way, one can establish uniform convergence along this sequence. 
To interpolate, one can employ the Avalanche Principle from [50] . 
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The paper [57j contains numerous applications of Theorem 4.10 to specific classes 



of subshifts, some of which will be described in Subsection |4.4| below. 

When [SSIISJ appeared, all known zero-measure results for Schrodinger operators 
defined by strictly ergodic subshifts (see, e.g., [H [H |251 [TOi [1101 EZl ) were 
covered by this approach, that is, they all held for subshifts that satisfy condition 
(B). Recently, Liu and Qu constructed examples of strictly ergodic subshifts that 
do not satisfy (B) but for which the associated Schrodinger operators do have 
zero-measure spectrum (and in fact the convergence of ^ log ||M£;.j^(r7,)|| to L{E) is 
uniform for every E G M.) fll6'. 

Naturally, once one knows that the spectrum has zero Lebesgue measure, one 
would like to determine its fractal (e.g., Hausdorff, lower and upper box counting) 
dimensions, as well as similar quantities such as thickness and denseness. These 
more delicate questions have been studied for a rather small number of examples, 
which will be discussed in subsequent sections. Zero-measure spectrum, on the 
other hand, is known in much greater generality, and this is the topic of the next 
subsection. 

4.4. Examples. In this subsection, we present several classes of popular examples 
of potentials that are ergodic, aperiodic, and take finitely many values (so that 
Kotani's central result applies) and discuss the validity of condition (B) (so that the 
associated Schrodinger operators have zero- measure spectrum). For more details, 
we refer the reader to [57] . 

4.4.1. Linearly Recurrent Subshifts and Subshifts Generated by Primitive Substitu- 
tions. A subshift ft over A is called linearly recurrent (or linearly repetitive) if there 
exists a constant K such that it v,w G yV(r2) with \w\ > K\v\, then is a subword 
of w. Clearly, every linearly recurrent subshift il satisfies (B). A popular way to 
generate linearly recurrent subshifts is via primitive substitutions. A substitution 
S : A ^ A* is called primitive if there exists fc £ N such that for every a,b G A, 
S'^{a) contains b. Such a substitution generates a subshift ft as follows. It is easy 
to see that there are m S N and a G A such that S"^{a) begins with a. If we iterate 
S"^ on the symbol a, we obtain a one-sided infinite limit, u, called a substitution 
sequence, fl then consists of all two-sided sequences for which all subwords are also 
subwords of u. One can verify that this construction is in fact independent of the 
choice of u, and hence is uniquely determined by S. Prominent example^ are 
given by 



at->- ab, b a 


Fibonacci 


ab, b ^ ba 


Thue-Morsc 


ai-^ ab, b ^ aa 


Period Doubling 


ai-^ ab, bh^ ac, c i-> db, d^ dc 


Rudin-Shapiro 



The following was shown in [72] (and independently in [66]): 

Proposition 4.11. // the subshift Q, is generated by a primitive substitution, then 
it is linearly recurrent and hence satisfies condition (B). 

It may happen that a non-primitive substitution generates a linearly recurrent 
subshift. An example is given by a i— >■ aaba, b ^ b. In fact, the class of linearly 



^These examples appear explicitly in many papers in the physics literature on Schrodinger 
operators generated by primitive substitution; compare [2l l6l 1141 1261 [78. .90..105.ril9| . 
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recurrent subshifts generated by substitutions was characterized in |58j In partic- 
ular, it turns out that a subshift generated by a substitution is hnearly recurrent if 
and only if it is minimal. 

4.4.2. Sturmian and Quasi- Sturmian Subshifts. Consider a minimal subshift ft over 
A. The (factor) complexity function p : Z+ — is defined by p{n) = #>Vn(n). 
Hedlund and Morse showed in [M] that O is aperiodic if and only if p{n) > n + 1 
for every n £ Z+ . Aperiodic minimal subshifts of minimal complexity, p{n) ~ n + 1 
for every n G N, exist and they are called Sturmian. If the complexity function 
satisfies p{n) — n + k for n > no, fc,no € N, the subshift is called quasi-Sturmian. 
It is known that quasi-Sturmian subshifts are exactly those subshifts that are a 
morphic image of a Sturmian subshift; compare 

There are a large number of equivalent characterizations of Sturmian subshifts; 
compare |17j . We are mainly interested in their geometric description in terms of 
an irrational rotation. Let a € (0, 1) be irrational and consider the rotation by a 
on the circle, Ra ■ [0, 1) [0, 1), Rad — {9 + a}, where {x} denotes the fractional 
part of X, {x} = X mod 1. The coding of the rotation Ra according to a partition 
of the circle into two half-open intervals of length a and 1 — a, respectively, is given 
by the sequences f„(a, 9) = X[o,q) (-^S^)- We obtain a subshift 

= {v{a,9) : e e [0,1)} = {v{a,9) : 9 e [0, 1)} U {w^^) (a) : fc G Z} C {0, if , 

which can be shown to be Sturmian. Here, t;i'^-*(a) — X(o,a] (-^S^'^O)- Conversely, 
every Sturmian subshift is essentially of this form, that is, if Q is minimal and has 
complexity function p{n) = n + 1, then, up to a one-to-one morphism, — Via for 
some irrational a € (0, 1). 

Using this description and some classical results in diophantine approximation, 
the following result was shown in |57j . 

Theorem 4.12. Every Sturmian subshift obeys the Boshernitzan condition (B). 

Moreover, establishing stability of (B) under morphic images, one obtains the 
following consequence. 

Corollary 4.13. Every quasi-Sturmian subshift obeys (B). 

4.4.3. Circle Maps. Let a € (0, 1) be irrational and /3 € (0, 1) arbitrary. The coding 
of the rotation R^ according to a partition into two half-open intervals of length 
/3 and 1 — /3, respectively, is given by the sequences w„(a,/3, 0) — X[o,;3) (^S^)- We 
obtain a subshift 

na,f} = {v{a,p,9):9^ [0,1)} C {0, 1}^. 

Subshifts generated this way are usually called circle map subshifts or subshifts 
generated by the coding of a rotation. 

The paper |57j established the following results for circle map subshifts in con- 
nection with property (B): 

Theorem 4.14. Let a € (0, 1) be irrational. Then the subshift Vta^p satisfies (B) 
for Lebesgue almost every /? e (0, 1). 



■^See also |69l 1701 11151 for results for Schrodinger operators arising from a specific class of 
non-primitive substitutions. 
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Theorem 4.15. Let a G (0, 1) be irrational with bounded continued fraction coef- 
ficients, that is, a„ < C. Then fla.p satisfies (B) for every (3 G (0, 1). 

Theorem 4.16. Let a e (0, 1) be irrational with unbounded continued fraction 
coefficients. Then there exists /3 G (0, 1) such that ^a.p does not satisfy (B). 

AAA. Interval Exchange Transformations. Subshifts generated by interval ex- 
change transformations (lETs) are natural generalizations of Sturmian subshifts. 
lETs are defined as follows. Given a probability vector A = (Ai,...,A„i) with 
\i > for 1 < i < m, we let /io = 0, /i^ = X]j=i'^j ' ^^'^ = [l^i-i, IJ-i)- Let 
T be a permutation of A^n = {1, • • • 7"^-}, that is, r G Sm, the symmetric group. 
Then A'^ = (At--i(i), . . . , \r-'^(m)) is also a probability vector, and we can form the 
corresponding jij and If . Denote the unit interval [0, 1) by /. The (A,t) interval 
exchange transformation is then defined by 

T : I ^ I, T{x) = X — fii^i + Mr(j)-i for x G /i, I <i < m. 

It exchanges the intervals li according to the permutation r. 

The transformation T is invertible, and its inverse is given by the (A'^,r~^) 
interval exchange transformation. 

The symbolic coding of x G / is LOnix) — i'\i T"{x) G li. This induces a subshift 
over the alphabet Am- ^x,t — {'-^i^) '■ x G /}. 

Sturmian subshifts correspond to the case of two intervals, as a first return map 
construction shows. 

Keane [100] proved that if the orbits of the discontinuities /^i of T are all infinite 
and pairwise distinct, then T is minimal. In this case, the coding is one-to-one and 
the subshift is minimal and aperiodic. This holds in particular if r is irreducible 
and A is irrational. Here, r is called irreducible if t({1, . . . , fc}) ^ ({1, . . . , k}) for 
every k < m and A is called irrational if the A^ are rationally independent. 

Regarding property (B), Boshernitzan has proved two results. First, in [20] the 
following is shown: 

Theorem 4.17. For every irreducible r G md for Lebesgue almost every \, 
the subshift J^a.t satisfies (B). 

In fact, Boshernitzan shows that for every irreducible r G Sm and for Lebesgue 
almost every A, the subshift VLx t satisfies a stronger condition where the sequence 
of n values for which r](n) is large cannot be too sparse. This condition is easily 
seen to imply (B), and hence the theorem above. 

In a different paper, [3T], Boshernitzan singles out an explicit class of subshifts 
arising from interval exchange transformations that satisfy (B) . The transformation 
T is said to be of (rational) rank k if the span a fc-dimensional space over Q (the 
field of rational numbers) . 

Theorem 4.18. IfT has rank 2, the subshift Q,\,r satisfies (B). 

4.5. Singular Continuous Spectrum. As seen in the previous section, the spec- 
trum has zero Lebesgue measure whenever condition (B) holds. This condition is 
satisfied by a wide class of models, in particular by all typical quasicrystal models. 
As pointed out in Theorem |4.10[ a consequence of zero-measure spectrum is the 
absence of absolutely continuous spectrum. That is, if <j{H^) has zero Lebesgue 
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measure, then (Ti^dHuj) ~ 0, since all spectral measures are supported by the spec- 
trum, and any measure supported by a set of zero Lebesgue measure must be purely 
singular by definition. 

To complement this, one can often show the absence of point spectrum. That 
is, there are a variety of tools that allow one to show that has no eigenvalues, 
and hence (Tpp(77^) = as well. Putting the two results together, one obtains that 
iJ^ has purely singular continuous spectrum. 

The primary tool that allows one to exclude eigenvalues is based on the Gordon 
lemma, which assumes that the potential has infinitely many suitably aligned local 
periodicities. Overall, this nicely implements the philosophy that aperiodic order is 
intermediate between periodic and random. The aperiodicity implies the absence 
of absolutely continuous spectrum via Kotani's theorem (and hence one does not 
have the spectral type that appears for a periodic medium) , while the order feature 
implies the absence of point spectrum via a fingerprint of local periodicity (and 
hence one does not have the spectral type that appears for a random medium). 

A potential y : Z — )■ M is called a Gordon potential if there are — > oo such 
that for every k, we have V{n) = V{n + qk) = V{n — qk) for 1 < n < q^. That is, V 
looks like a periodic potential around the origin, as one sees at least three suitably 
aligned periodic unit cells there, and the period may be chosen arbitrarily large. 
The following Gordon lemma is based in spirit on 82J. In this particular form it 
was shown in |67) . 

Lemma 4.19. Suppose V is a Gordon potential. Then, for every E, the difference 
equation 

u{n + 1) + u{n - 1) + V{n)u{n) = Eu{n) 

has no non-trivial square- summable solutions. In particular, the associated 
Schrddinger operator H in ^^(Z), given by 

[Htl;]{n) ^ il}{n + 1) + i>{n - 1) + V{n)iP{n), 

has empty point spectrum. 

By ergodicity, T-invariance, and the Gordon lemma, if one can show that 

/Lt ({a; € 57 : VLj is a Gordon potential}) > 0, 

then 

/i ({cj G ^l : has empty point spectrum}) = 1. 

On the other hand, for any aperiodic minimal subshift, at least one of its elements 
fails to have the required Gordon three- block symmetries [ID]. Thus, one cannot use 
this appraoch to show uniform absence of eigenvalues for all a; € fl. Nevertheless, 
results to this effect are known, established with the following variant of the Gordon 
lemma. 

Lemma 4.20. Suppose : Z — > M is such that there are q^ ^ oo such that 
V(n) — V{n + gfe) for I < n < qi^. Suppose further that E is such that 



(4) sup 

k 



< oo. 



Then, the difference equation 

u{n + 1) + u{n - 1) + V{n)u{n) = Eu{n) 
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has no non-trivial solutions that are square- summable on Z+ and hence E is not an 
eigenvalue of the associated Schrddinger operator H in ^^(Z). In particular, if the 
assumption Q holds for every E in the spectrum of H , then H has empty point 
spectrum. 

In many quasicrystal models, the existence of hierarchical structures gives rise 
to a so-called trace map, which in turn can often be used to ensure that Q holds 
for all energies in the spectrum. Thus, the analysis then reduces to finding suitable 
squares of arbitrary length starting at the origin. 

Thus, in the symbolic setting at hand, the observations above give rise to prob- 
lems that concern the subword structure of the potentials, and hence fall in the 
general area of combinatorics on words. 

Let us describe the results that have been obtained in this way for the examples 
discussed above. In all these results, the choice of the sampling function / is more 
restricted than above. Namely, one usually assumes that /(w) — g{u!o) with an 
injective map g ; — > M. 

4.5.1. Subshifts Generated by Primitive Substitutions. Suppose S* is a primitive sub- 
stitution over the alphabet A and let C be the subshift associated with it. 
Recall that it is strictly ergodic and denote the unique invariant probability mea- 
sure by /i. The index of f2 is given by the largest fractional power occurring in 
an (and hence any) element of r2. Formally, the index is defined as follows. Given 
w S Wo and any prefix v of w, we denote the word w'^v with k G Z+ by w^, 
where r = fc -f y^. Then, indf;(u') — sup{r G Q n [l,oo) : £ Wo} and 

ind(17) = sup{indo(w) : w € Wo}- 

The following result was shown in [iO] 

Theorem 4.21. If S is a primitive substitution and the associated subshift Q sat- 
isfies ind(ri) > 3, then fj, ({cj G i7 : is a Gordon potential}) > 0. Consequently, 
H {{ui G 17 : has empty point spectrum}) ~ 1. 

The underlying idea is simple. Since the subshift is invariant under S, any word 
appearing with index strictly greater than 3 generates by iteration of S" a sequence 
of words whose lengths go to infinity and whose index is bounded away from 3. 
This allows one to bound from below the frequency with which third powers occur 
and hence yields measure estimates on the Gordon three-block conditions that are 
good enough to show that the limsup of these sets must have positive measure. 
Since the elements of the limsup of these sets give rise to Gordon potentials, the 
result follows. 

Applications of this theorem include the Fibonacci substitution (since ind(r2) > 
indo(a6aa6) > 3.2), the period doubling substitution (since ind(r2) > indo(a6) > 
3.5), and many others. Of course, the result does not apply to the Thue-Morse 
substitution, which is famous mainly because ind(r2) — 2. Unfortunately, it is still 
open whether the point spectrum is almost surely empty in the Thue-Morse case. 
The Gordon approach fails due to small index, and no other methods are known 
that yield an almost sure result 

^See [39| for a precursor dealing with the period doubling case. In this special case it was later 
shown that the absence of eigenvalues even holds for all a; £ O I42J . 

^The absence of eigenvalues for a dense set of oj's can be established in this example and 
many others using palindromes instead of powers f89| . However, using palindromes one cannot 
prove the absence of eigenvalues for a full measure set of a;'s |66| . 
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4.5.2. Sturmian and Quasi- Sturmian Subshifts. Damanik, Killip, and Lenz showed 
the foUowing result in |52j (see also |53| for a uniform result for almost every Stur- 
mian subshift). 

Theorem 4.22. For every Sturmian subshift fJ, H^^ has empty point spectrum for 
every w G fJ. 

This result was the culmination of a sequence of partial results. Among those, 
we single out Siito [145) . who proved empty point spectrum for one a and one to, 
Bellissard et al. |13|, who proved it for all a and one cj, Delyon-Petritis [6^, who 
proved it for almost all a and almost all w, and Kaminaga |98| . who proved it for 
all a and almost all w. Here, a € (0, 1) \ Q denotes the slope associated with a 
Sturmian subshift. Recall that Sturmian subshifts are in one-to-one correspondence 



with (0, 1) \ Q. Here, [H [51 HI Ei] used Lemma |4:20| whereas ^ ^ used 
Lemma 14.191 

The extension of Theorem |4.22| to the quasi-Sturmian case was obtained by 
Damanik and Lenz in |55| . 

Theorem 4.23. For every quasi-Sturmian subshift f2, has empty point spec- 
trum for every w G fi. 

4.5.3. Circle Maps. Recall that a circle map subshift is determined by the param- 
eters as (0, 1) \ Q and /3 € (0, 1). It is strictly ergodic and we denote the unique 
ergodic measure by /i. Delyon and Petritis proved the following in |67) . 

Theorem 4.24. For almost every a and every j3, the corresponding circle map 
subshift D, is such that fJ.{{uj G : 14; is a Gordon potential}) = 1. Consequently, 
fi {{uj £ Q : has empty point spectrum}) = 1. 

In fact, the full measure set of a values is explicitly described in terms of the 
continued fraction expansion. The condition was weakened by Kaminaga in [98], 
still however excluding an explicit zero measure set. This weaker condition was only 
shown to imply /i {{u G 17 : T4; is a Gordon potential}) > 0, which of course is still 
sufficient to allow one to deduce ^ {{lo G : H^^ has empty point spectrum}) = 1. 

4.5.4. Interval Exchange Transformations. Recall that an lET subshift is deter- 
mined by an irreducible permutation r and a probability vector A. Cobo et al. 
showed the following result in [33] (see also [68|). 

Theorem 4.25. For every irreducible permutation t and almost every \, the as- 
sociated lET subshift n is such that /i {{io G 12 : is a Cordon potential}) = 1. 
Consequently, /i ({w G 11 : H^^ has empty point spectrum}) — 1. 

4.6. Transport Properties. Quasicrystal models have behavior that is markedly 
different from the periodic and random cases in many different respects. In the 
previous subsections we have seen that the spectrum is typically a zero-measure 
Cantor set, while for periodic and random potentials it is always given by a finite 
union of non-degenerate compact intervals. Moreover, the spectral type is typically 
singular continuous, while it is always absolutely continuous in the periodic case 
and pure point in the (one-dimensional) random case. 

In this subsection we consider yet another perspective from which the quasicrys- 
tal model behavior is expected to differ from the behavior of the periodic and 
random cases. Namely, we consider the spreading of wave packets under the time- 
dependent Schrodinger equation. That is, given a Schrodinger operator and a 
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normalized element ^p of ^^(Z), we consider = e~**"^'^-0, where e~^*^'-' is de- 
fined via the spectral theorem. Then, V'(') satisfies the time-dependent Schrodinger 
equation idt^j{t) = Hi^iplt) with initial condition ip{0) = 4'- The quantum me- 
chanical interpretation is that the probability of finding the quantum particle at 
site n e Z at time t e M is given by a{n,t) := \{Sn,tp(t))\'^ ■ The initial state is 
naturally localized in some fixed compact set, up to a small error, since it belongs 
to £'^{Z). More specifically, one is often interested in the initial state ip — Sq (or 
some Sn), which is completely localized. After having fixed the initial state, one 
is then interested in how fast ipit) spreads out in space, or more specifically, how 
long one has to wait until a{n,t) is no longer negligibly small at some distant site 
n. In general, this is a difficult problem. Questions of this kind are easier to study 
for compound quantities; that is, some averaging in n and/or t helps one generate 
quantities for which interesting statements can be proven. 

A popular way to average in time is to consider Cesaro averages, 

a{n,T)^^[ a{n,t)dt=^[ | (<5„, ^i- 

Let also 



Mp{t) = J2i^ + \nnain,t), Mp(T) ^(1 + |nna(n, T), p > 0. 

Notice that for t (resp., T) fixed, a{-,t) and a{-,T) are probability distributions on 
Z, and hence the quantities above are (1 plus) the p-th moment of the respective 
probability distribution. Here we assume that the initial state is either a Dirac 
delta function or at least sufficiently well localized so that these moments are finite. 

Wave packet spreading then is reflected by growth in time of these moments. To 
detect power-law growth, one introduces the so-called transport exponents 

= hmsup ^°SM,{t) ^ ^ ^ ^.^.^^ logM.jt) ^ 

t^oo Plogt <^oo plogt 

T->oo plogT T-)-oo plogT 

Each of thsee four functions of p is non-decreasing in p and takes values in the 
interval [0,1]. 

In view of the monotonicity of the transport exponents, it is natural to consider 
their limiting values for small and large values of p. Thus, denote 

af = lim/3''=(p), a J = lim/?='=(p), 

af = \imP'^{p), a J = lim P'^ip). 

We note that there are other useful ways of capturing wave packet spreading, and 
refer the reader to |65j for a comprehensive survey. 

The transport exponents take the constant value for random potentials and 
(at least the time-averaged quantities) the constant value 1 for periodic potentials. 
Thus if one is able to prove the occurrence of fractional values of the transport 
exponents, one exhibits wave packet spreading that is strictly intermediate between 
the periodic and random cases. Results of this kind are notoriously difficult to 
establish. The few known results for quasicrystal models will be described in detail 
in later sections on Fibonacci and Sturmian potentials. 
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5. The Fibonacci Hamiltonian 

The Fibonacci Hamiltonian is the most prominent model in the study of elec- 
tronic properties of quasicrystals. It is given by the discrete one-dimensional 
Schrodinger operator 

(5) [Hx,uju]{n) — u{n + 1) -|- u(ri — 1) + Ax[i-a,i) ("-a -f co mod l)u(n), 

where A > is the coupling constant, a — is the frequency, and uj € [0, 1) 

is the phase. An alternative way to obtain the same potential is via the Fibonacci 
substitution; see Section |4.4.1| above. This operator family has been studied in 
many papers since the early 1980's (see, e.g., pi Ffgl 1501 lOTl 15^ ITM [TUil [TII51 [1151 
1126111271 1150] for early works in the physics literature), and numerous fundamental 
results are known. In this section we describe the current "state of the art" for this 
model. 

5.1. Trace Map Formalism. Even the earliest papers on the Fibonacci Hamil- 
tonian realized the importance of a certain renormalization procedure in its study, 
see [1031 1126j . This led in particular to the consideration of a certain dynamical sys- 
tem, the so-called trace map, whose properties are closely related to many spectral 
properties of the operator ([s]). The existence of the trace map and its connection 
to spectral properties is a consequence of the invariance of the potential under a 
substitution rule. This works in great generality; see [3l I41j and references therein. 

The one-step transfer matrices associated with the difference equation Hx^^u — 
Eu are given by 



Tx^u{m,E) 



~ -^Xfi-Q.i) ("^ct + ™od 1) — 1 
1 



Denote the Fibonacci numbers by {-Ffc}, that is, Fq = Fi = 1 and F^+i — Fk + Fk-i 
for k > 1. Then the fact that the potential for zero phase is invariant under the 
Fibonacci substitution implies that the matrices 

M_i(i?)=(^J -^^y Mo(i?)=(^f 

and 

MkiE) = Tx,o{Fk,E) X • • • X r;,,o(l,^) for fc > 1 
obey the recursive relations 

Mk+AE) ^ Mk-iiE)MkiE) 
for fc > 0. Passing to the variables 

XkiE) = ^TrMfe(S), 

this in turn implies 

(6) Xk+i{E) = 2xk{E)xk-i{E) - Xk-2[E) 

for fc > 1, with x^i{E) = 1, xo(-E) = E/2, and xi ^ {E - A)/2. The recursion 
relation ([6]) exhibits a conserved quantity; namely, we have 

(7) Xk+i{Ef + Xk{Ef + Xk-i{Ef - 2xk+i{E)xk{E)xk-i{E) - 1 = ^ 
for every fc > 0. 
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Given these observations, it is then convenient to introduce the trace map 

(8) T -.R^ ^ M^ T{x, y, z) = {2xy - z, x, y). 

Aside from the context described here, this map appears in a natural way in 
problems related to dynamics of mapping classes |79], Fuchsian groups [18], 
number theory [19], Painleve sixth equations [29l |95], the Ising model for qua- 
sicrystals [21 [551 11531 1154] . the Fibonacci quantum walk |137[ 1138) . and others 
[3 [Sni 11481 1155j . See [55] or [H] for an algebraic explanation of this universality. 
We refer the reader also to |139l I140j for further reading on the Fibonacci trace 
map. 

The following function, 

G{x, y, z) — x'^ + y'^ + z"^ — 2xyz — 1, 

is invariant under the action of Tj^ and hence T preserves the family of cubic 
surfacef[3 

(9) Sx = |(a;, y, z) eR^ : x^ + y^ + z^ - 2xyz = 1 + ^| . 

Plots of the surfaces 6*0.01 and 6*0.5 are given in Figures [ij and [2j respectively. 
Denote by £\ the line 

It is easy to check that £\ C 5a. 

Siito proved the following central result in [145'. 

Theorem 5.1. An energy E belongs to the spectrum of Hx^uj */ o,i^d, only if the 
positive semiorbit of the point ( ^2 ^ ' T 1 ^^der iterates of the trace map T is 
bounded. 



'The function G{x,y,z) is usually called the Fricke character or Fricke-Vogt invariant. 
The surface So is called the Cayley cubic. 
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To obtain this theorem, Siito argued as follows. Denote 

ak = {EeR:\xk{E)\<l} 

and 

These sets depend on the coupling constant A, and whenever we want to make 
this dependence explicit, we will write cTk^x and Sfc.A- An analysis of the trace 
recursion (|6| shows that the sets are decreasing, and hence it is natural to 
consider their limit S = H ^fe- Clearly, if G S, then {xniE)} remains bounded 
due to ([7]). On the other hand, the analysis of the trace recursion ^ also yields that 
whenever E ^ Tik for some fc, then \xn-kiE) \ obeys an explicit super-exponentially 
growing lower bound. That is, the sequence {xn{E)} remains bounded if and only 
a E ef:. Notice that the point {^^, f , 1) is just {xi{E),xoiE),x^i{E)), so that 



Theorem |5.1| follows as soon as S = S is established. The inclusion ECS holds 
since a^. is precisely the spectrum of the canonical periodic approximant of period 
Fk and the fact that these periodic approximants converge strongly. The inclusion 
E C I] holds since one can use the boundedness of {xn{E)} for e S along with 
the Gordon lemma to show that no solution for this energy is square-summable at 
+CXD, which implies that E must be in the spectrum. 



5.2. Hyperbolicity of the Trace Map. Let f : M ^ M be a diffeomorphism 
of a Riemannian manifold M. Let us recall that an invariant closed set A of the 
diffeomorphism / is hyperbolic if there exists a splitting of a tangent space T^M = 
E^ © i?" at every point x € A such that this splitting is invariant under Df, and 
the differential Df exponentially contracts vectors from stable subspaces {E^} and 
exponentially expands vectors from unstable subspaces {i?"}. A hyperbolic set A 
of a diffeomorphism / : M — >■ M is locally maximal if there exists a neighborhood 
[/(A) such that 

A = fl /"([/). 

We will consider diffeomorphisms of a surface, dim M = 2, and hyperbolic sets of 
topological dimension zero. In this case a locally maximal hyperbolic set A can be 
locally represented as a product of "stable" and "unstable" Cantor sets and C". 
Both Cantor sets C and are dynamically defined in the sense of the following 
definition. 

Definition 5.2. Let / C be a closed interval. A Cantor set C <Z I is dynam- 
ically defined if there are strictly monotone contracting maps ipi,^2, ■ ■ ■ li^k '■ I ^ 
1,^,(1) nip J (I) ^$ ifiy^j, such thatC = n„eN-f", where h = V'lW U • • • U 
and In+i = U • • • U ipk{In)- 

If "01, V'2i ■ • ■ ; "0*; are C^'''^-functions, then the Cantor set has zero measure, de- 
pends continuously on tpi^ ... ^tjjk, and is "regular" in many other ways. We will be 
interested in the Hausdorff dimension and the thickness of the Cantor sets C" and 
C" . Denote the Hausdorff dimension of the set C by dimjj C. 

In our case, dim// A = dim// C + dimj/ C"; see [1221 [ISQ]. Moreover, if / e C"^ 
depends smoothly on a parameter, then dim// A is also a smooth function of the 
parameter; see [120) . 
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Definition 5.3. Let C C M now be an arbitrary Cantor set and denote by I its 
convex hull. Any connected component of I\C is called a gap ofC. A presentation 
of C is given by an ordering hi ~ {Un}n>i of the gaps of C. If u d C is a 
boundary point of a gap U of C , we denote by K the connected component of 
I\{Ui U C/2 U • • • U Un) (with n chosen so that U„ = U) that contains u and write 

With this notation, the thickness t(C) and the denseness 9{C) of C are given 

by 

(10) r(C) = supinf r(C,Z^,M), e{C) = inf supr(C,Z^, w). 

Tlie tliiclcness and tlie denseness of a Cantor set C are related to tiie Hausdorff 
dimension of C by tlie following inequalities (of. [1291 Section 4.2]), 

(11) , < dimn C < 
log(2+^)- « -log(2+5(i^) 

For more details on thickness, see [7TJ 11231 I129j . An important property of 
thickness was discovered by Newhouse |125j : 

Theorem 5.4. If Ci and C2 are two Cantor sets and t(Ci) • t(C2) > 1, then the 
sum Ci + C2 contains an interval. In the special case Ci = C2 ='■ C , we have that 
t{C) > 1 implies that C + C is an interval. 

Consider the restriction Tx : 5a — ?> Sx of the trace map T from ^ to the 
invariant surface Sx, Tx — T\s^. Denote by U,x the set of points in Sx whose full 
orbits under Tx are bounded. 

Theorem 5.5. For every A > 0, the set ^x is a locally maximal hyperbolic set of 
Tx '. Sx ^ Sx- It is homeomorphic to a Cantor set. 



Theorem 5.5 was proved for A > 16 by Casdagli [3Tj, for small values of A by 
Damanik and Gorodetski |37], and finally for all A > by Cantat [55]. 

It is easy to check that ^x C Sx, and therefore the set of points on Ix whose 
forward semiorbits are bounded is exactly equal to lx^W''{Q.x). Then the spectrum 
Yjx is affine equivalent to the set Ix H W^{Vtx). It is known that £x intersects the 
leaves of W''{^x) transversally for A > 16 (see Section 2 in [5^), and for sufficiently 
small A > (see [17]). This allows one to consider the spectrum as a dynamically 
defined Cantor set for these values of the coupling constant. Therefore the following 
holds. 

Corollary 5.6. For A > 16, as well as for sufficiently small A > 0, the spectrum 
Hx is a dynamically defined Cantor set, and hence: 

(i) For every small £ > and every x G cr{Hx.u}), we have 

AuTiH {{x ~ e, X + e) n a{Hx,u:)) = dim^ {{x - e,x + e)n cf{IIx^ui)) 

= dimn <j{IIx,uj) 
= dims a{IIx,uj)- 

(ii) The Hausdorff dimension dmiH <y{IIx,uj) is an analytic function of X, and 
is strictly between zero and one. 

We conjecture that the same holds for all A > 0. 
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5.3. Quantitative Characteristics of the Spectrum at Large Coupling. 

The fact that the box counting dimension of the spectrum exists and coincides 
with its HausdorfF dimension allows one to determine the asymptotic behavior of 
this A-dependent quantity in the large coupling limit. In fact, Damanik-Embree- 
Gorodetski-Tcheremchantsev proved the following in [46] . 

Theorem 5.7. We have 

lim (diniEA) • log A = log(l + V2). 

A— f oo 

Let us briefly explain how this result is obtained. Recall that the spectrum is 
related to the spectra of the canonical periodic approximants by 



Sa = Pi Sfc_A = Pi CTfc.A U cTfc+i.A- 



fc>i fe>i 

Since each periodic spectrum ak,x is a finite union of non-degenerate compact in- 
tervals and the lengths of these intervals can be shown to be decaying, it is natural 
to use Sfc,A as one possible cover of T,\ and estimate the HausdorfF dimension of 
Sa from above in this way. On the other hand, since each interval of ak,\ can be 
shown to have non-empty intersection with E^, one can estimate the box counting 
dimension of T,\ from below in this way. We observe how crucial it is that these 
dimensions coincide here. Thus, the analysis of the participating intervals comes 
down to proving good estimates for their length. 

To estimate the length, one makes use of the following basic fact from one- 
dimensional Floquet theory. The preimage of the open interval (—1,1) under xj- 
consists of exactly F^. disjoint open intervals, on which Xk is strictly monotone. In 
fact, in this particular case, the same statement is true for the corresponding closed 
intervals (i.e., the periodic spectra in question have all their gaps open). Thus, the 
length of one of these intervals (say / — [a, b]) can be estimated as follows. Since 

2 = \xk{a)-xk{b)\ = I \x'^{E)\dE, 



we have 



< 1/1 < 



maxsg/ \x'^(E)\ min^g/ \x'^.{E)\' 

In order to prove estimates for |a::'^(£')|, one differentiates the trace recursion (|6| 
and proceeds inductively, making use of the trace invariant (j?]). This approach was 
pioneered by Raymond |136j and then used in many subsequent papers. In this 
inductive approach, it turns out to be important to determine, for a given energy 
E in one of the intervals of ak,\, in how many of the earlier sets (Jk',\i k' < fc, 
the energy E in question lies. This gives rise to a combinatorial question that was 
completely answered in [46] . Combining these combinatorial results with the length 
estimates one can prove in this way for the intervals in question, the overall strategy 
above yields the following specific estimates: 

log(l + V2) 



(12) dimn Sa < 



log(i [(A-4) + V(A-4)2-12]) 



for A > 8, 



(13) dim^ Ea > for A > 4. 

^ ' ^ log (2A + 22) 
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Theorem |5.7| is then a direct consequence of these estimates and the fact that the 
Hausdorff dimension and the box counting dimension of S]> are equal for A > 16. 

5.4. Quantitative Characteristics of the Spectrum at Small Coupling. 

Fractal properties of Ea for small A were studied in '^SJ. Among many other 
things, that paper established the following pair of theorems. 

Theorem 5.8. We have 

lim dim T,\ = 1. 

\->-o 

More precisely, there are constants Ci, C2 > such that 

1 - CiA < dimEA < 1 - C2A 
for A > sufficiently small. 
Theorem 5.9. We have 

lim t(T,x) — 00. 
More precisely, there are constants C3 , C4 > such that 

C3A-1 < r(EA) < 0(Ea) < C4A-1 
for A > sufficiently small. 



Theorem 5.8 is a consequence of the connection (11) between the Hausdorff 
dimension of a Cantor set and its denseness and thickness, along with the estimates 
for the latter quantities provided by Theorem |5.9| 

Let us briefly explain how Theorem |5.9| can be obtained. The surface Sq, 
from ([9]), is called the Cayley cubic; it has four conic singularities and can be 
represented as a union of a two dimensional sphere (with four conic singularities) 
and four unbounded components. The restriction of the trace map to the sphere 
is a pseudo-Anosov map (a factor of a hyperbolic map of a two-torus), and its 
Markov partition can be presented explicitly (see [3T] or [371 US])- For small values 
of A, the map T : S\ S\ "inherits" the hyperbolicity of this pseudo-Anosov map 
everywhere away from the singularities. The dynamics near the singularities must 
be considered separately. Consider the dynamics of T near one of the singulari- 
ties, say, near the point p = (1, 1, 1). The set Per2(T) of periodic orbits of period 
two is a smooth curve that contains the point p and intersects S\ at two points 
(denote them by pi(A) and P2(A)) for A > 0. Finite pieces of stable and unstable 
manifolds of pi{X) and P2(A) are a distance of order A from each other. In order to 
estimate the thickness (and the denseness) of the spectrum Sa, we notice first that 
the Markov partition for T : Sq Sq can be continuously extended to a Markov 
partition for T : S\ ^ Sx. The extended Markov partition is formed by finite parts 
of the stable and unstable manifolds of pi(A), P2(A), and the other six periodic 
points that are continuations of the three remaining singularities. Therefore the 
size of the elements of these Markov partitions remains bounded, and the size of 
the distance between them is of order A. The natural approach now is to use the 
distortion property (see, e.g., [129] ) to show that for the iterated Markov partition, 
the ratio of the distance between the elements to the size of an element is of the 
same order. The main technical problem here is again the dynamics of the trace 
map near the singularities, since the curvature of S\ is very large there for small 
A. Nevertheless, one can still estimate the distortion that is obtained during a 
transition through a neighborhood of a singularity and prove boundedness of the 
distortion for arbitrarily large iterates of the trace map. This implies Theorem|5.9| 
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5.5. The Density of States Measure. Let us now turn to the formulation of 
results involving the integrated density of states, a quantity of fundamental impor- 
tance associated with an ergodic family of Schrodinger operators. The integrated 
density of states (IDS) was introduced in Section 3.2 in a more general context, 
and represents the distribution function of a density of states measure - a measure 
supported on the spectrum and, in particular, reflecting the asymptotic distribution 
of eigenvalues of finite dimensional approximations. 

Denote the density of states measure of the Fibonacci Hamiltonian for a given 



coupling constant A by dN\. Repeating the definition from Section 3.2 in this 
particular case, we have 

, , , , #|eigenvalues of Hx ^ ri ,,1 that are < E\ 
(14) Nx{E) = lim ■ ■ ' = — -, 

where i?A.t^,[i,n] is the restriction of H\^i^ to the interval [1, n] with Dirichlet bound- 
ary conditions, and the limit does not actually depend on the phase w. 

It is interesting to analyze the regularity of the density of states measure. This 
question was studied for general potentials [33 1371 1771 1114' , random potentials 
[23 [Hg, and analytic quasi-periodic potentials [3 [22l |23l [24, ^ ,811 |83l [MIJ . In 
the case of Fibonacci Hamiltonian, the IDS is Holder continuous. 

Theorem 5.10. For every A > 0, there exist C\ < 00 and > such that 

\NxiEi) - Nx{E2)\ < Cx\Ei ~ E^^' 

for every Ei, E2 with \Ei — E2\ < 1. 

This follows directly from |52]; see also [38l [54l [92l [93l [96] for some previous 
related results. 

It is also interesting to obtain the asymptotics of the optimal Holder exponent 
for large and small couplings. 

In the large coupling regime, we have the following [51] : 

Theorem 5.11. (a) Suppose A > 4. Then for every 

31og(a-i) 
^ 21og(2A-|-22)' 

there is some S > such that the IDS associated with the family of Fibonacci 
Hamiltonians satisfies 

\Nx{Ei) - Nx{E2)\ < \E1-E2r 
for every Ei, E2 with \Ei — E2\ < 5. 
(b) Suppose A > 8. Then for every 

1 > 



31og(a-i) 



12 



21og(i((A-4) + /(A^2 

and every < S < 1, there are Ei,E2 with < \Ei — £^2! < S such that 
\Nx{Ei) - Nx{E2)\ > \E1-E2r. 

Thus the optimal Holder exponent is asymptotically ^'2 fog a ^ large cou- 

pling regime. 

The proof is based on the self-similarity of the spectrum and an analysis of the 
periodic approximants (in the spirit of the proof of Theorem 5.7). 
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In the small coupling regime, we have the following |51j : 

Theorem 5.12. The integrated density of states N\{-) is Holder continuous with 
Holder exponent where 7a — > | as \^ Q, and 7a < | for small A > 0. 
More precisely: 

(a) For any 7 G (0, \), there exists Aq > such that for any A G (0, Aq), there 
exists (5 > such that 

\Nx(Ei) - Nx{E2)\ < \E1-E2r 

for every Ei, E2 with \Ei — E2\ < S; 

(b) For any sufficiently small A > 0, there exists 7 = 7(A) < ^ such that for 
every d > 0, there are Ei, E2 with < \Ei — -E2I < <^ (^nd 

\N^{E^) ~ N^{E2)\ > 1^1 -SsT- 

The proof uses the trace map formalism and a relation between the IDS of Hx i^ 
and the measure of maximal entropy for the trace map Tx. Namely, the density 
of states measure is proportional to the projection (along the stable manifolds) to 
£\ of the normalized restriction of the measure of maximal entropy fJ-ms^x{T\) to 
an element of the Markov partition. After that, the proof uses a comparison of 
expansion rates of Tx and Tq (and is reminiscent of the proof of Holder continuity 
of conjugacies between two hyperbolic dynamical systems). 

Another interesting feature of the Fibonacci Hamiltonian is the uniform scaling 
of the density of states measure. Namely, the following holds [50l: 

Theorem 5.13. There exists < Ao < c» such that for A G (0, Aq), there is 
dx G (0, 1) so that the density of states measure dNx is of exact dimension dx, that 
is, for dNx-almost every £' G M, we have 

lim^°g^-(^-^'^ + ^^^dA. 

e^O log £ 

Moreover, in (0, Aq), c^a is a function of X, and 

limc^A = 1- 

The proof is based on the relation between dNx and fJ.ma.x{Tx), and the exact 
dimensionality of hyperbolic measures |H1 1108L 1133) . 

The Hausdorff dimension of the spectrum is an upper bound for dx, but a priori 
it is not clear whether these numbers must coincide. Barry Simon conjectured that 
for a large class of models these quantities must be different. The next result shows 
that this conjecture is true for the weakly coupled Fibonacci Hamiltonian [50] . 

Theorem 5.14. For A > sufficiently small, we have dx < dim/f Sa- 

The proof is based on the comparison of the measure of maximal entropy for Tx 
(which is "responsible" for dx) and the equilibrium measure for the potential given 
by minus the log of the expansion rate. The Hausdorff dimension of the unstable 
projection of the latter is equal to dim// Sa, and the thermodynamical description 
of this measure (see |122j ) implies that for any other ergodic invariant measure, the 
dimension of its unstable projection is strictly smaller. In order to prove that those 
two measures are actually different, one uses the fact that the measure of maximal 
entropy is an equilibrium measure that corresponds to zero potential. Therefore it is 
enough to show that the two potentials under consideration are not cohomological. 
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which can be done using a comparison of multipliers of different periodic orbits 
of Ta. 

5.6. Gap Opening and Gap Labeling. The spectrum jumps from being an 
interval for A = to being a zero-measure Cantor set for A > 0. Hence, as the 
potential is turned on, a dense set of gaps opens immediately. It is natural to ask 
about the size of these gaps; see [Hj. These gap openings were studied in [S] for 
the Thue-Morse potential (where the gaps open as a power of A) and in [11] for 
the period doubling potential (where some gaps open linearly, and some others are 
superexponcntially small in A). In the Fibonacci case, all gaps open linearly |49| : 

Theorem 5.15. For A > sufficiently small, the boundary points of a gap in the 

spectrum T^x depend smoothly on the coupling constant A. Moreover, given any 
one-parameter continuous family {[/a}a>o of gaps of "Ex, we have that 

y \UX\ 

exists and belongs to {0,oo). 

Theorem 1 5 . 1 5| follows again from dynamical properties of the trace map. Namely, 
each singularity of the Cayley cubic S'o gives birth to two periodic points on the 
surface ^a, A > 0. The distance between the periodic points is of order A. The 
stable manifolds of these periodic points "cut" gaps in £x that correspond to gaps in 
the spectrum. The curves formed by the families of the periodic points are normally 
hyperbolic manifolds of the trace map, and hence (see [551 1135j ) their strong stable 
manifolds form a foliation. This implies that the size of each gap is also of order 



A (as A — >■ 0), and Theorem 5.15 follows. 

The limit in Theorem 5.15 certainly depends on the family of gaps chosen. In 
order to study this dependence, one needs to use some labeling of the gaps. As is 
well known, the density of states produces such a gap labeling. That is, one can 
identify a canonical set of gap labels, which is only associated with the underlying 
dynamics (in this case, an irrational rotation of the circle or the shift-transformation 
on a substitution-generated subshift over two symbols) , in such a way that the value 
of N{E, A) for i? € K\Sa must belong to this canonical set. In the Fibonacci case, 
this set is well-known (see, e.g., [12l Eq. (6.7)]) and the general gap labeling theorem 
specializes to the following statement: 

(15) {N{E, X) : E eR\T.x} C {{ma} : to € Z} U {1} 

for every A ^ 0. Here {ma} denotes the fractional part of ma, that is, {ma} = 
ma — [tockJ . Notice that the set of gap labels is indeed A-independent and only 
depends on the value of a from the underlying circle rotation. Since a is irrational, 
the set of gap labels is dense. In general, a dense set of gap labels is indicative 
of a Cantor spectrum and hence a common (and attractive) stronger version of 
proving Cantor spectrum is to show that the operator "has all its gaps open." For 
example, the Ten Martini Problem for the almost Mathieu operator is to show 
Cantor spectrum, while the Dry Ten Martini Problem is to show that all labels 
correspond to gaps in the spectrum. The former problem has been completely 
solved [4], while the latter has not yet been completely settled. Indeed, it is in 
general a hard problem to show that all labels given by the gap labeling theorem 
correspond to gaps, and there are only few results of this kind. It turns out that the 
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stronger (or "dry" ) form of Cantor spectrum holds for the weakly coupled Fibonacci 
Hamiltonian [15] : 

Theorem 5.16. There is Aq > such that for every A € (0, Aq], all gaps allowed 
by the gap labeling theorem are open. That is, 

(16) {N{E, X) : E eR\T.x} ^ {{ma} : m e Z} U {1}. 

Complete gap labeling for the strongly coupled Fibonacci Hamiltonian was shown 



by Raymond in |136j . where he proves ( 16 1 for A > 4. We conjecture that ( 16 1 holds 
for every A > 0. 

Using the gap labeling, we can refine the statement of Theorem |5.15| For m e 
Z\ {0}, denote by Um{^) the gap of Ea where the integrated density of states takes 
the value {ma}. Then, the following result from [49] holds: 

Theorem 5.17. There is a finite constant C* such that for every m G Z \ {0}, 

\UmiX)\ Cm 

hm — r— — = - — - 
A^o |A| \m\ 

for a suitable Cm G (0, C*). 



To see why Theorem 5.17 holds, notice that each family of gaps converges (as 
A — > 0) to a point of intersection of io with a stable manifold of one of the singu- 
larities. The intersections that have larger labels are in a sense "produced" from 
intersections with smaller labels by the action o f the inverse of the trace map. For 
gaps with small labels, we know from Theorem |5.15| that HmA_^o ^^^xp^ < C* for 



some constant C* > 0. The length (in coordinates on the two-torus covering Sq) 
of the piece of the stable manifold from the singularity to the point of intersection 

after k applications of the map is of order ^ ^^2^ ) ^ 1"^!' ^^'^ contraction 
that will be applied to the gap is of order 



V5-l\ /\/5-lV-(^) 1 



5.7. Square and Cubic Fibonacci Hamiltonians. Since spectral questions for 
Schrodiger operators in two (and higher) dimensions are hard to study, it is natural 
to consider a model where known one-dimensional results can be used. In particular, 
let us consider the Schrodinger operator 

(2) 

[H)^ il;]{m, n) =%l){m -f 1, n) -f 1/1(771 — 1, n) -f -(/(m, n -|- 1) -|- ip{m, n — 1)+ 
+ A {x[i~a.i){'ma mod 1) + X[i-a.i){na mod 1)) -0(m, n) 

in ^^(Z^). The theory of tensor products of Hilbert spaces and operators then 
implies that o'{H^ ) = + Sa- This operator and its spectrum have been studied 
numerically and heuristically by Even-Dar Mandel and Lifshitz in a series of papers 
[751 [Til 175] (a similar model was studied by Sire in [144) 1. Their study suggested 
that at small coupling, the spectrum of Ea + Ea is not a Cantor set; quite on the 
contrary, it has no gaps at all. 

It turns out that this is indeed the case [15]: 

(2) 

Theorem 5.18. For A > sufficiently small, cr{H\^ ) = Ea + Ea is an interval. 
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The same statement holds for the cubic Fibonacci Hamiltonian (i.e., the anal- 
ogously defined Schrodinger operator in P'il?) with spectrum + + 
Theorem |5 . 18| follows from the estimates for the thickness of Ea from Theorem |5.9| 
and Newhouse's Gap Lemma (Theorem 5.4). Section 7.3 shows numerical illustra- 
tions of the finite approximations S^^a + '^k,\ and S^^a + '^k,\ + '^k.\^ along with 
an exploration of the number of disjoint intervals that make up these sets. 



5.8. Transport Properties. There is a substantial number of papers that inves- 
tigate the transport exponents associated with the Fibonacci Hamiltonian; see, for 
example, [ISllMlllSllMllSllSnilSIliniElIMlIMl IM] ■ While we won't describe 
all the known results, we want to at least highlight some of them and put them in 
perspective. As pointed out earlier, one of the fascinating features of quasicrystal 
models is that their aperiodic order being intermediate between periodic and ran- 
dom is reflected in a number of ways, be it in a spectral way (by spectral measures 
being purely singular continuous) or in a transport behavior way. Here we want 
to address the latter point. All the papers listed above have to goal of proving 
estimates that show that the transport properties of the Fibonacci Hamiltonian are 
markedly different from those of periodic or random media. 

Since there is ballistic transport (all transport exponents are equal to one) in 
the periodic case and no transport (all transport exponents are equal to zero) in 
the random case, one therefore wants to show that the transport exponents take 
values in the open interval (0, 1). Proving non-trivial lower bounds turns out to be 
comparatively easier and was accomplished in the late 1990's [3H1IM] for zero phase. 
Several subsequent papers then went on to extend the lower bound to all phases 
and improved the estimates [351 ISH EHl ED IS21 [S3] ■ Upper bounds for transport 
exponents, on the other hand, proved to be elusive for some time. Note a key 
difference here: to bound transport exponents from below, one "only" has to show 
that some portion of the wave packet moves sufficiently fast. On the other hand, 
to bound transport exponents from above, one essentially has to control the entire 
wave packet and show that it does not move too fast (i.e., ballistically) . Thus, it is 
potentially easier to prove upper bounds on transport that are dual to the type of 
lower bound that had been established, and this indeed turned out to be the case. 
The papers [331 llOlj showed that at least some non-trivial portion of the wave 
packet moves slowly. Full control and hence genuine upper bounds for transport 
exponents were finally obtained in 2007 and later [HI [63l E] . 

Let us now state some of the transport results explicitly. Some general remarks 
that should be made are the following: 

(a) Almost all results conc ern t ime-averaged quantities (i.e., the exponents 
/J*(p) defined in Section 4.6). 



(b) Most papers focus on the case ?/'(0) — (5o. We will limit our attention here 
to this case as well. 

(c) The optimality of the known estimates improves when p and/or A are large. 
In particular, the bounds are known to be tight in the limit A,p f oo. 

(d) For finite values of A and p, the method of choice to obtain the best known 
bound varies. 

(e) For A and p large enough, the transport exponent may exceed the dimension 
of the spectrum. 
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Here is a result from [62] that establishes the best known estimates for zero phase 
and given A and p: 

Theorem 5.19. Suppose A > and set 



7 = Dlog(2 + V8 + A2) 
{where D is some universal constant) and 

''"^°^[201og(l + «) 

Then, the time- averaged transport exponent corresponding to the initial state "0(0) 
Sq and zero-phase Fibonacci Hamiltonian i?A,o obey 

P+2k 



(17) 



13^ {p) > 



(p+i)(7+K+i/2)' p<27 + i; 

.4t' p>27 + 1. 



Here is a result from [M| that concerns the regime of large A and p: 
Theorem 5.20. Consider the Fibonacci Hamiltonian Hx u cii^-d the initial state 



V'(O) = 6o. For A > v24, we have 



.± 21og(l + a) 



and for A > 8, we have 



log(2A + 22) ' 
21og(l + a) 



12 



log(i [(A-4) + y(A^ 
Both estimates holds uniformly in uj. In particular, 

hm • log A = 2 log(l + a), 

A— J-oo 

and convergence is uniform in u). 

In fact, the upper bound can be proved also for the non-time-averaged quantities, 
as shown in |64j . 

Theorem 5.21. Consider the Fibonacci Hamiltonian H\^^ and the initial state 
'0(0) = 5q. For A > 8 and uniformly in lo, we have 

^± < 21og(l + a) 

" log (i [(A - 4) + V(A-4)2-12]) ■ 

6. Sturmian Potentials 

The Fibonacci potential is a special case of a Sturmian potential. The latter are 
obtained if a in the definition of the potential, V{n) — Ax[i-q,i) (f^ct -I- wmod 1), is 
a general irrational number in (0, 1). The Fibonacci case corresponds to the choice 

Given an irrational a e (0, 1), consider its continued fraction expansion 

1 



1 



0.1 



02 



as 
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with uniquely determined G 1^+- Truncating the continued fraction expansion 
of a after k steps yields the rational number Pk/lki which is the best rational 
approximant of a with denominator bounded by Qk+i — 1- The following recursions 
hold: 

Pk+i = ak+iPk +Pk-i. po = 0, pi = l, 

^fc+i = ak+iQk + qk-i- qo = 1, 9i fli- 

(In the Fibonacci case a = , we have = 1 a,ndpk/qk = Fk-i/Fk.) A number 
of the results for the Fibonacci Hamiltonian described in the previous section have 
been generalized to the Sturmian case under suitable assumptions on the continued 
fraction coefficients {a^}. In this section, we explain what these results are, and 
how the proofs had to be modified. 

6.1. Extension of the Trace Map Formalism. Let us the denote the discrete 
Schrodinger operator on i'^(Z) with potential V{n) = Ax[i_q i) (rta + wmod 1) by 
Hx^a,u- Strong approximation again shows that the spectrum of H\ a,ui does not 
depend on w, and may therefore be denoted by SA,a- The one-step transfer matrices 
associated with the difference equation = Eu are given by 

^^^^ ^(^^^^^ j^£;-Ax[w)(-" + -nrodl) -l^j ^ 
The matrices 

and 

Mk{E)^Tx,c..o{qk.E) x ••• x Tx,cA^,E) for A: > 1 
obey the recursive relations 

Mk+i{E) - Mk-i{E)Mk{EY''+' 
for fc > 0; see [T31 Proposition 1]. Passing to the variables 

Xk{E) = ^TrAffe(S), 

this in turn implies via the Cay ley-Hamilton theorem that Xk+i{E) can be ex- 
pressed as an explicit function of (suitable Chebyshev polynomials applied to) 
Xk{E),Xk-i{E),Xk-2{E) for A: > 1; see [13l Proposition 2]. These recursion re- 
lations exhibit the same conserved quantity as before; namely, with 

Xk+iiE) = ^TiiMk{E)Mk-i(E)), 

we have 

ik+i{Ef + Xk{Ef + Xk-i{Ef - 2S:k+i{E)xk{E)xk-i{E) - 1 = ^ 
for every fc > 0; see [131 Proposition 3]. 
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6.2. Results Obtained via an Analysis of the Trace Recursions. Notice that 
the key difference with the Fibonacci case is that, in general, the sequence of traces 
may not be obtained by iterating a single map. In this sense, there is in general 
no direct analog of the trace map. However, as we have just seen, the underlying 
structure of recursive relations extends nicely. The substitute for the dynamical 
analysis of the Fibonacci trace map will have to lie in studying the dynamics of an 
initial point under the successive application of a sequence of maps, the elements of 
which are dictated by the continued fraction expansion of a. These developments 
are still in their early stages. In the following we will concentrate on the known 
results that can be established by simply exploiting the recursive relations, without 
employing sophisticated tools from dynamical systems theory. 

The first result that establishes a clean analogy with the Fibonacci case is the 



following analog of Theorem 5.1 which was established in [inj 



Theorem 6.1. Fix A > and a e (0, 1) irrational. An energy E belongs to the 
spectrum Ea.q if and only if the sequence {xk{E)} is bounded. 



The proof of Theorem |6.1| follows the same line of reasoning as the proof of 
Theorem |5.1[ which was outlined in the previous section. In particular, one obtains 
canonical covers of the spectrum, which are useful in the estimation of its dimension. 
Let us make this explicit. As before, define the sets 

'^x,c.,k^{EeR:\xkiE)\<l} 

and 

The same reasoning shows that the sets SA,a,fc are decreasing in k and the spectrum 
is the limiting set, that is, 

k>l 

see [ISl Proposition 4]. 



A refinement of this description of the spectrum in the Sturmian case due to 
Raymond 1 136] allowed Liu and Wen to obtain the following estimates for the 
Hausdorff dimension of the spectrum in the large coupling regime |118) . 

Theorem 6.2. Suppose A > 20 and a G (0, 1) is irrational with continued fraction 
coefficients {ak}. Denote 

M = liminf(ai • • •flfc)^/'' G [l,oo]. 

/c— >-oo 

(a) If M = oo, then dim// S^.q = 1. 

(b) If M < oo, then dim// „ belongs to the open interval (0, 1) and obeys the 
estimates 

2 log A/ + log3 
21ogAf-log-^ 

log M - log 3 1 
101og2 - 31og ' logM - log y2(^ J ■ 
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A study of the box counting dimension of SA,a was carried out in the follow-up 
paper |7B] by Fan, Liu, and Wen. Among other things, they showed that for A > 20, 
the Hausdorff dimension and the box counting dimension of Sa.q coincide whenever 
the sequence {a^} is eventually periodic. 

The transport exponents in the Sturmian case were studied in the papers f5Sl 
[SH [501 EH I121j . The following result from ^62] gives dynamical lower bounds for 
all values of A and p, provided a has bounded continued fraction coefRcients. 

Theorem 6.3. Suppose A > and a G (0, 1) is irrational with ak < C. With 
7 = D log(2 + + A2) • lim sup - ^ 



k=l 

{where D is some universal constant) and 

^ logjVU/A) 

the transport exponents associated with the operator H\,a.o and the initial state 
1^(0) = ^0 obey 

P+2k 



p<2a + l; 
p>2a + l. 



The following result from [121] gives dynamical upper bounds in the large cou- 
pling regime. 

Theorem 6.4. Suppose A > 20 and a G (0, 1) is irrational with continued fraction 
coefficients {ak\ and corresponding rational approximants {pk/lk}- Denote 

D = lim sup loggfc. 

fc— i-oo 

Then, the transport exponents associated with the operator H\ afi and the initial 
state "0(0) — 6q obey 

+ 2D 



log^ 

Moreover, if ak > 2 for all k, then 

~+ D 



7. Numerical Results and Computational Issues 

In this section, we provide numerical illustrations of a number of the results 
described in this survey. These calculations focus on the Fibonacci Hamiltonian, 
though many could readily be adapted to the Sturmian potentials described in 
the last section. We begin by studying approximations to the spectrum for the 
Fibonacci model in one dimension, then investigate estimates of the integrated 
density of states based on spectra of finite sections of the operator. Finally, we 
address upper bounds on the spectrum in two and three dimensions. In all cases, 
we set the phase w to zero. 
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7.1. Spectral Approximations for the Fibonacci Hamiltonian. We begin by 
calculating the spectrum for the kth periodic approximations to the Fibonacci 
potential. The analysis described in Section [5] suggests several ways to compute (7^, 
which turn out to have varying degrees of utility. 

Given a candidate energy E, one can test if _E e ct^ by iterating the trace 
recurrence ^ and testing if |a;fe(i5)| < 1. In principle, this simple approach enables 
investigation for arbitrarily large values of k. However, two key obstacles restrict 
the utility of this method: (i) it does not readily yield the entire set Ufe; (ii) as k 
increases, the intervals that comprise ak become exponentially narrow, beyond the 
resolution of the standard floating point number system in which such calculations 
are typically performed. However, this approach can yield some useful results, 
particularly in the small coupling regime where the decay of the interval widths is 
most gradual, or when one is only interested in some narrow set of energy values. 
(This method of calculation was used to produce illustrations in [35].) 

To obtain the entire set , one might instead use the recurrence ^ to construct 
the degree--F'fc polynomial XkiE), then determine the regions where |xfe(i?)| < 1 by 
finding the zeros of the polynomials Xk{E) + 1 and Xk{E) — 1 using a standard root- 
finding algorithm. For all but the smallest k this approach is untenable. Coefficients 
of Xk{E) grow exponentially in k; e.g., for A = 4, 

xeiE) = i£;i3 - 16E^^ + '^E^^ - 1616£;i" + l^^E^ - 16272£;« + ISSSOS^ 

+ 20160£;s - 37133£;5 - 17056^;'* + ^^E^ + 25104£;2 + ^^E + 560. 

The magnitude of these coefficients, compounded by the proximity of the roots for 
larger values of A and fc, leads to inaccurate root calculations, a phenomenon well 
studied by numerical analysts; see, e.g. |1241I151) . Indeed, it is not uncommon for 
the computed roots to be so inaccurate as to have significant spurious imaginary 
parts. 

There is a more robust approach to computing the approximate Fibonacci spec- 
trum ak ■ One can view ak as the exact spectrum of a related Schrodinger operator 
with a potential having period F^. The spectrum of this operator is the union of 
Fk non-degenerate intervals whose end points are given by the eigenvalues of the 
two i^fc -dimensional matrices Jk+ and Jfe_: 

/ t'l.fe 1 ±1 \ 



Jk± — 



with unspecified entries beyond the tridiagonal section equal to zero; see, e.g., [1491 
Ch. 7]. Here the potential values Vn,k are given by 

(18) Vn,k = Ax[i-_Ffc_i/Ffe,i)("-F'fe-i/-F'fe mod 1). 

This approach, which we use for the computations described below, has also been 
employed in the context of Fibonacci computations by Even-Dar Mandel and Lif- 
shitz [73] . 

The standard procedure for computing all the eigenvalues of a symmetric ma- 
trix begins by applying a unitary similarity transformation to reduce the matrix to 
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symmetric tridiagonal formj^ Floating-point arithmetic introduces errors into this 
process, resulting in the exact tridiagonal reduction of a matrix that differs from 
the intended matrix by a factor that scales with the precision of the floating point 
arithmetic system, the coupling constant A, and the dimension Fk. The eigenvalues 
of this tridiagonal matrix are then approximated to high accuracy via a procedure 
known as QR iteration }131j . Remarkably, this iteration does not introduce sig- 
nificant errors beyond those incurred by the reduction to tridiagonal form; for a 
discussion of this accuracy, see pi ll52| . Overall, this process requires 0{F^) floating 
point arithmetic operations and the storage of 0{F^) floating point numbers. (The 
conventional procedure for reducing the matrix to tridiagonal form destroys the 
zero structure present in Jk±-) Of course, the upper estimate T,k.x — crk.x U cr^+i^x 
then requires computation of all eigenvalues of four matrices. 

Significant insight can be gleaned from numerical calculations involving small 
to moderate values of k. For example. Figure |3] shows a^.x for A G [0, 2] and 
k = 1, . . . , 8, while Figure [4] shows the upper bounds T,k,x for the same range of A 
and k. Since Sg^A = os.x U ag^x, for A > the spectrum is the union of 34 and 55 
intervals. 

To develop conjectures (e.g., regarding dim one would like use approxima- 
tions to Yix for larger k. Two fundamental challenges arise: (i) the 0{F^) work and 
0{F^) storage becomes prohibitive; (ii) while non-degenerate, the intervals in ak,x 
become exponentially small and exponentially close together. This phenomenon 
is illustrated in Figure [5] The utility of the numerical results degrade when the 
size of these bands and gaps approaches the order of the error in the numerical 
computation|^ On contemporary commodity computers, computations of S^^a up 
to roughly k — 20 (requiring all eigenvalues of matrices of dimension F20 = 10,946 
and F21 = 17,711) is feasible, provided A is sufficiently small for the results to 
be accurate. Recently Puelz has proposed an improved approach that ameliorates 
challenge (i) above by reducing the required work to 0(F|) and storage to 0{Fk), 
and challenge (ii) by using extended precision arithmetic |134) . 

To estimate the box-counting dimension of T,x (assuming it exists), we use the 
definition 

A- fc\ V logCs(e) 

dmiB(A) = hm — — — , 

£^0 logl/e 

where Cs(e) counts the number of intervals of width e that intersect S, 

Csie) #{j e Z : [je, (j + l)e) n 5 ^ 0}. 

Note that dmiB{Y,k.x) — 1 for all fc, since Sfc^A is the union of finitely many closed 
intervals. Still, one gains insight into dimB(EA) from log(Cs^ ^ (£))/log(l/£) for 
finite values of e and various fc, as can be seen in Figure |6] For fixed A, the 
resulting estimates of dimBCSx) (taken, e.g., as infEg(o,i) log(Ci;^ (e))/ log(l/£)) 
apparently improve as k increases; lower values of k are suitable for larger values of 
A. However, with this approach it is difhcult to accurately estimate the critical value 



^Methods such as the Lanczos algorithm excel at computing a few eigenvalues of large sym- 
metric matrices 11311 Ch. 13]. These methods are not feasible here, for all eigenvalues of Jj,^ are 
required. However, if one is only interested in a narrow band of energies, these methods can be 
highly effective. 

® More subtly, the formula \l8\ incurs significant rounding errors for large n and k, resulting 
in errors on the diagonal of Jf^± of size A. For greater accuracy, one should use the equivalent 
formulation v„ k = -'*X[Ffc -Fi._i .Fj.) ("--'^fe— 1 rnod ^fe), which is more robust. 
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at which diniB(I]A) = 1/2. (A rough estimate, suggested from Figure [6j is A ~ 4; 



see the discussion preceding Problem 8.7 below.) More accurate approximations 



will require computations with larger values of k than are feasible with the method 
described above. 

Finally, Figure [7] explores numerical computations of the thickness, defined 
in (10). As established in Theorem 5.9 the thickness t(Ea) behaves like 1/A as 
A 4, 0. As A decreases we see this behavior mirrored in the upper bounds S/c^a, up 
to some point where T(T,k.\) rapidly increases: Sfc.A is the union of no more than 
Fk + Fk+i intervals separated by gaps that diminish as A J, 0. 



7.2. Density of States for the Fibonacci Model. We next turn to an inves- 
tigation of the exponent of Holder continuity of the integrated density of states 
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Figure 4. The upper bound Y^k.x = (^k,x U (Jk+i.\ on the Fi- 
bonacci spectrum, as a function of A £ [0, 2]. For k = 8 and A = 2, 
Sfe,A is the union of 42 disjoint intervals. 



(IDS) for the Fibonacci model, discussed in Section 5.5 To estimate Nx{E) in 



equation (14), one must compute all the eigenvalues of -ffA.[i,n]j the restriction of 
H\ to sites [1, n] with Dirichlet boundary conditions. This restriction is an n x n 
tridiagonal matrix; because this matrix lacks the corner entries present in Jk± in 
the last subsection, its eigenvalues can readily be efficiently computed for large val- 
ues of n (say n < 10^ on contemporary desktop computers). While computational 
complexity is no longer such a constraint, accuracy still is: for large n and A, some 
eigenvalues of i?A,[i,ri] are closer than the precision of the floating point arithmetic. 
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Figure 5. Exponential decay of the largest intervals and small- 
est gaps in the approximations Sfc.A, for coupling constants A — 
2, 4, 8, 16, 32, as computed in MATLAB's double-precision floating 
point arithmetic. The data points that are plotted in gray are 
likely dominated by computational errors. 






Figure 6. Estimates of AimB{^\) for various values of A, based 
on the upper bounds S^^a for various k. The dashed horizontal line 
denotes log(l -I- A/2)/log(A), to which dmiBi'i^x) tends as A — ?> oo 
( Theorem |5.7[ ). The gray horizontal lines in the bottom plots show 
the upper and lower bounds (12)-(13|. 
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Figure 7. Thickness of S^.a as a function of A for three values 
of k, consistent with Theorem 15.91 



rendering, for example, \Ei — = for theoretically distinct eigenvalues Ei and 
E2 of i/A,[l,„][] 

Figure |8] shows estimates of the IDS based on computations with n — 10,000 for 
A values ranging from the trivial case of no coupling (A = 0) to strong coupling 
(A = 8). The fine structure of the spectrum is evident in Figure|9j which repeatedly 
zooms in upon subsets of the spectrum of the finite section -ffA,[i,Ti] for A = 1 and 
n = 100,000. (The numerical concerns discussed in the last paragraph do not affect 
these figures.) 

We now explore the Holder continuity of the integrated density of states. In 



consideration of (14), define 



Nn,xiE)= Hm 



^{eigenvalues of ^^A.[i.ri] that are < E} 



Figure 10 investigates the large A behavior of the Holder exponent addressed in 
Theorem 5.111 based on computations with finite sections of dimension n = 10,000. 



Indeed, we see asymptotic behavior like '^^2 fog a ^ ' ^^'^ moreover the figure suggests 
that the dimension of the measure is smooth in this regime. 

7.3. Spectral Estimates for Square and Cubic Fibonacci Hamiltonians. 

As described in Section 5.7 the estimates Sfc,A for the one-dimensional Fibonacci 



spectrum can readily be translated into approximations for the square and cubic 
cases, as investigated by Even-Dar Mandel and Lifshitz [73]. As described in The- 
orem |5.18[ Sa need not be a Cantor set, especially for small coupling constants. 
This behavior is apparent in Figures 11 and 12 which illustrate Efc,A + ^k.x and 



Sfc^A + Sfe^A + Sfc. A for various values of k and A. For a finite range of small A values, 
the spectra are intervals that branch into a greater number of intervals as k and A 
increase. Figure [13] shows the growth in the number of intervals present in these 
approximations as a function of A for three different values of k. This plot makes 
evident rapid (but not always monotone) growth in the number of intervals with 



^'^By its structure, must have n distinct eigenvalues. However, scenarios where eigen- 

values are exceptionally close are well-known in the numerical analysis community; see, e.g., 
Wilkinson's W:^^ matrix [ml Sec. 7.7]. 
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Figure 8. Approximations to the integrated density of states for 
the Fibonacci model with various values of the coupling constant, 
A, based on n = 10,000. 



A. Figure [M] illustrates the opening and closing of gaps for the square problem, 
revealing an intriguing structure for finite k. How does this structure develop as k 
increases, and, indeed, is it reflected in Ea + At present these questions remain 
open. 

Tables[l]and[2]investigate the square and cubic spectral estimates more precisely, 
giving the values of A where multiple intervals first emerge. These results confirm 
and sharpen the observation of Even-Dar Mandel and Lifshitz [73] that T,k,\ + Efe,A 
transitions from one to two intervals near A = 1.3, while S^^^ + ^k,\ + 5]fe,A makes 
the same transition near A = 2. For these finite values of fc, it is apparent that 
Efc^A + ^k,\ and S^^^ + ^k,x + ^k,x both transition to two intervals, then three 
intervals, and so on. What do these calculations suggest about the limit k — )• oo? 
For example, is the A value at which S^^a + '^k,\ transitions from two to three 
intervals converging? Is there a finite span of A values for which Yi^^x + Tik^x persists 
as the union of two intervals as fc — >■ oo, or does S^ + E^ transition from one interval 
directly to a Cantorval or Cantor set? (See Problems 8.7 and 8.8 below.) 
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Figure 9. Eigenvalues of Hx.[i^n] for A = 1 and n = 100,000, 
drawn as vertical lines to aid visibility. The first plot shows the 
entire spectrum; the gray boxes denote the region on which the 
next plot zooms. 
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Figure 10. Illustration of Theorem 5.11 based on numerically- 
computed eigenvalues from finite sections -ffA,[i,n] for n = 10,000. 
Here 5 — 0.025 and the minimization is over Ei,E2 G a{Hx\i^n])- 



Let Xk^m denote the value of A at which Tj^ x + '^k.x (or Sj. ^ + Ej. ^ + Ti^ x) 
first splits from m to m + 1 intervals as A increases, with Afe^o — 0. (Our detailed 
computations suggest that, for small values of m, there is only one such point of 
transition; for larger numbers of intervals, gap closings complicate the picture, as 



seen in Figure 14 ) Figure 15 plots \k,m. — Afc_m-i as a function of fc for m = 
1, . . . , 7 for the square and cubic Hamiltonians. Do the transition points converge 
as fc — ^ oo? First consider the plot on the left, for the square Hamiltonian. For 
m = 1 and m = 2, \k^m appears to converge; however, the points of transition to 
TO > 3 intervals do not show such consistency: it is unclear if these \k,m values 
are converging. It may be that the coupling constants at which T,x,k + '^\,k breaks 
into TO > 3 intervals are converging to the point at which the spectrum breaks into 
TO = 3 intervals as fc — ?> oo. Now consider the plot on the right of Figure [15) for the 
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Figure 11. Approximations S^.a + '^k.\ to the spectrum of the 
two-dimensional Fibonacci operator, as a function of A. For fc = 8 
and A = 4, I]fc_;)v-|-Efe_;vis the union of 311 disjoint intervals. 

cubic Hamiltonian. In contrast to the square case, these results suggest the \k.ra 
values converge to distinct points as fc — ^ cx) for all values m = 1, 2, . . . , 7 shown, 
inviting the conjecture that there exist A values for which + + is the union 
of TO disjoint intervals for all to > 1. 



8. Conjectures and Open Problems 

In this final section we discuss various open problems that are suggested by the 
existing results and address generalizations, strengthenings, and related issues. 

We begin with open problems for the Fibonacci Hamiltonian. The existing quan- 
titative results concern estimates for dimensional properties of the spectrum, the 
density of states measure, and the spectral measures, as well as estimates for the 
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Figure 12. Approximations E^.a + ^k.\ + ^k.\ to the spectrum 
of the three-dimensional Fibonacci operator, as a function of A. 
For k = 7 and A = 7, T,k,\ + ^k,\ + ^k,\ is the union of 482 disjoint 
intervals. 



transport exponents. In almost all cases, the asymptotic behavior is known in the 
regimes of small and large coupling. While the bounds we obtain are monotone, 
we would like to understand whether the quantities themselves have this property: 

Problem 8.1. Are the various quantities we consider [in particular, dim// T,x) 
monotone in X? 

The known estimates for dimensional properties of the spectral measures are 
the only ones that are clearly not optimal, and in particular do not identify the 
asymptotics in the extremal coupling regimes. For the density of states measure, 
which is an average of spectral measures, we have much better information. Can 
one find ways to find equally good estimates for spectral measures? 
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A A 
Figure 13. Number of intervals in the spectral approximations 
5^fe,A + '^k,\ and Sfc,A + Sfe,A + 5]fc,A, as a function of A. 

Table 1. Estimates of A*, the A value for which the thickness 
of Sfc^A equals one, along with Xk,rm the coupling constant where 
Sfe_A + ^k,x splits from m to m + 1 intervals, for m = 1, . . . , 4. 



k 






Afe,2 


Afe,3 


•^fe,4 


6 


1 


313172936 


1.313172936 


1 


624865906 


1.649775155 


1 


708521471 


7 


1 


298964798 


1.298964798 


1 


543759898 


1.548912772 


1 


596682038 


8 


1 


296218739 


1.296218739 


1 


494856217 


1.514291562 


1 


520122025 


9 


1 


294303086 


1.294303086 


1 


445808095 


1.492410878 


1 


512965310 


10 


1 


293935333 


1.293935333 


1 


442778219 


1.446787662 


1 


472813609 


11 


1 


293679331 


1.293679331 


1 


430901095 


1.436192692 


1 


437915282 


12 


1 


293630242 


1.293630242 


1 


402035016 


1.415460742 


1 


426586813 


13 


1 


290031553 


1.293596081 


1 


392730451 


1.412863780 


1 


419815054 


14 


1 


288819456 


1.293589532 


1 


382510414 


1.404399139 


1 


408704405 


15 


1 


287431935 


1.293584975 


1 


380466052 


1.399646887 


1 


400190389 


16 


1 


287269802 


1.293584102 


1 


380121550 


1.388518687 


1 


397593470 


17 


1 


287084388 


1.293583494 


1 


379851608 


1.387310733 


1 


395556145 


18 


1 


287062735 


1.293583377 


1 


379806139 


1.385835331 


1 


393702258 


19 


1 


287037977 














20 


1 


287035086 















Problem 8.2. What can one say about the spectral measures? In particular, are 

their dimensional properties uniform across the hull and/or across the spectrum? 
Moreover, what are the asymptotics as A ^ and A f oo ? 

We know that dimS^ goes to one as A goes to zero. In addition, we would be 
interested in the following: 
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Figure 14. Approximations S^^a + '^k,\ to the two-dimensional 
Fibonacci spectrum Yi\ + Ti\ as in Figure 11, magnified to show 



the opening and closing of gaps as A increases. How this structure 
affects + Ti\ is not presently understood. 

Table 2. Estimates of Xk,m, the coupling constant where S^^a + 
Sfe A + A splits from m to m + 1 intervals, for to = 1, . . . , 4. 



k 


Afe,i 




Aa;,3 


Aa;,4 


6 


2 


025741216 


2 


544063632 


2 


573539294 


2 


842670115 


7 


2 


012664501 


2 


438240772 


2 


511570744 


2 


606841186 


8 


2 


011113604 


2 


376933028 


2 


498126298 


2 


498926850 


9 


2 


009524869 


2 


364541039 


2 


435665993 


2 


473875055 


10 


2 


009337409 


2 


357357667 


2 


412613336 


2 


421115367 


11 


2 


009145619 


2 


355932060 


2 


399696274 


2 


408616763 


12 


2 


009123008 


2 


355107791 


2 


392573154 


2 


401253561 


13 


2 


009099880 


2 


354944739 


2 


391094663 


2 


397036745 


14 


2 


009097154 


2 


354850520 


2 


390282080 


2 


393347062 


15 


2 


009094365 


2 


354831891 


2 


390113912 


2 


393329303 


16 


2 


009094036 


2 


354821128 


2 


390021550 


2 


393302392 


17 


2 


009093700 


2 


354819000 


2 


390002443 


2 


393300376 




Problem 8.3. Does the right- derivative oj dm\Y,\ exist at zero? 

On a technical level, an answer to the following question is crucial: 

Problem 8.4. Show that l\ is transversal to W^{Q\) for all X values. 

A proof of this property would immediately result in an extension of several 
results that are known to hold for small (and/or large) values of A to all values of 
A. 

Let us now turn to the higher-dimensional separable analogs of the Fibonacci 
Hamiltonian (e.g., the square or cubic Fibonacci Hamiltonian) . Recall that the 
spectrum of such an operator is given by the sum of the one-dimensional spectra, 
which in turn are Cantor sets. Recall also that at sufficiently small coupling, these 
sum sets are intervals, while at sufficiently large coupling, they are Cantor sets as 
well. Concretely, this uses that if the thickness of a Cantor set C is larger than 1, 
then C -|- C is an interval by Theorem |5.4| and, on the other hand, if the upper box 
counting dimension of C is strictly less than 1/2, then C -I- C is a Cantor set. It 
is natural to ask what shape the higher-dimensional spectra have at intermediate 
coupling, that is, we wish to study how the transition from C + C being an interval 
to being a Cantor set happens when the thickness of C decreases. 

Definition 8.5. A compact set C C is a Cantorval if it has a dense interior 
{i.e., int{C) = C), it has a continuum of connected components, and none of them 
is isolated. 

Here is a general result on the occurrence of Cantorvals in the context of taking 
sums of Cantor sets |123) : 

Theorem 8.6. There is an open set hi in the space of dynamically defined Cantor 
sets such that for generic Ci, C2 G U, the sum Ci + C2 is a Cantorval. 

Unfortunately, this result does not provide any specific and verifiable genericity 
conditions that would allow one to check that the sum of two given specific Cantor 
sets is indeed a Cantorval. Thus, for our purpose we need a solution to the following 
problem. 
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Problem 8.7. Provide specific verifiable conditions on a Cantor set C which imply 
that the sum C + C is a Cantorval. 

Ideally, such a criterion would be applicable to the spectrum of the Fibonacci 
Hamiltonian and establish that, say, the spectrum of the square Fibonacci Hamil- 
tonian is a Cantorval for intermediate values of the coupling constant A. The next 
step would then be to study the transitions between the three regimes. We ask 
whether there are two sharp transitions; compare [24 for closely related numerical 
evidence and discussion. 

(2) 

Problem 8.8. Let H)^ be the separable square Fibonacci Hamiltonian. Prove that 
there are values < A' < A" < oo such that for A G (0, A'), the spectrum a{H\) is 
an interval {or a finite union of intervals) , for X Cz {X' , X"), it is a Cantorval, and 
for X G (A", oo), it is a Cantor set. 

Notice that this will provide an example of a (topologically!) new structure of 
the spectrum for "natural" potentials. 

Moving on from the Fibonacci case, which has a description via a substitution 
rule as well as via a simple quasi-periodic expression, there are two natural choices 
of a more general setting. 

For a different choice of the underlying substitution rule, one always has an as- 
sociated trace map. However, our understanding of the dynamics of such a trace 
map is in general far more limited than the one in the Fibonacci case. As a conse- 
quence, outside of the Fibonacci case there is a scarcity of quantitative results for 
dimensional issues (such as the dimension of the spectrum, the dimension of the 
density of states measure, or the dimension of the spectral measures). For example, 
here is a simple open problem that is currently completely out of reach: 

Problem 8.9. Study other trace maps {e.g., period doubling and Thue-Morse); in 
particular, find the asymptotics of the HausdorfJ dimension of the spectrum as the 
coupling constant tends to zero or infinity. 

The other natural generalization of the Fibonacci potential is to replace the 
golden ratio in its quasi-periodic description by a general irrational number. Thus, 
we discuss some open problems for Sturmian potentials next. 

Let us say that two Cantor sets Ci and C2 on are diffeomorphic if there 
are neighborhoods [/i(Ci), 1^2(^2), and a diffeomorphism f ■ Ui ^ U2 such that 
/(Ci) = C2. 

Problem 8.10. Suppose that a — [ai, 02, . . .] and /3 = 62, . . .] are such that for 
some fc e Z and all large enough i € Z+ we have bi^j^ = a^. Prove that in this case, 
the Sturmian spectra Sa.q Cii^'d '^\,I3 diffeomorphic. 

Notice also that due to the ergodicity of the Gauss map, a solution of this problem 
would also imply that the following long standing conjecture is correct: 

Problem 8.11. For any fixed A > 0, the dimension dim// is almost every- 
where constant in a. 

Finally, let us emphasize that most of the questions related to higher dimensional 
models (described in Section [S]) are completely open. So we formulate an extremely 
general problem: 
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Problem 8.12. Study spectral properties {e.g., the shape of the spectrum, and the 
type of the spectral measures) and transport properties of higher dimensional opera- 
tors; for example, study these questions for the particular case of the Laplacian on 
the graph associated with a Penrose tiling. 
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